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We present medium-resolution (A/AA ~ 20, 000) ultraviolet spectra covering the 1 155-1760 A spectral range of the Seyfert 1 galaxy 
Mrk 509 obtained using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST). Our observations were 
obtained simultaneously with a Low Energy Transmission Grating Spectrometer observation using the Chandra X-ray Observatory, 
and they are part of a multiwavelength campaign in September through December 2009 which also included observations with 
XMM-Newton, Swift, and INTEGRAL. Our spectra are the highest signal-to-noise observations to date of the intrinsic absorption 
components seen in numerous prior ultraviolet observations. To take advantage of the high signal-to-noise ratio, we describe special 
calibrations for wavelength, flat-field and line-spread function corrections that we applied to the COS data. We detect additional 
\& . complexity in the absorption troughs compared to prior observations made with the Space Telescope Imaging Spectrograph (STIS) on 

HST. We attribute the UV absorption to a variety of sources in Mrk 509, including an outflow from the active nucleus, the interstellar 
medium and halo of the host galaxy, and possible infalling clouds or stripped gaseous material from a merger that are illuminated by 
the ionizing radiation of the active nucleus. Variability between the STIS and COS observation of the most blue-shifted component 
(#1) allows us to set an upper limit on its distance of < 250 pc. Similarly, variability of component 6 between FUSE observations 
limits its distance to < 1.5 kpc. The absorption lines in all components only partially cover the emission from the active nucleus with 
covering fractions that are lower than those seen in the prior STIS observations and are comparable to those seen in spectra from the 
Far Ultraviolet Spectroscopic Explorer (FUSE). Given the larger apertures of COS and FUSE compared to STIS, we favor scattered 
light from an extended region near the active nucleus as the explanation for the partial covering. As observed in prior X-ray and UV 
spectra, the UV absorption has velocities comparable to the X-ray absorption, but the bulk of the ultraviolet absorption is in a lower 
ionization state with lower total column density than the gas responsible for the X-ray absorption. We conclude that the outflow from 
the active nucleus is a multiphase wind. 
<& '. 
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1 . Introduction ment of the intergalactic medium (IGM). The most distant AGN 

known, quasars at redshifts of 6 or more, have abundances 



Outflows from galaxies powered by an active galactic nucleus at sol ar or super-sola r levels (lHamann & Ferlan J fl993l [l99 



(AGN) may play an important role in the chemical enrich- iPentericci et all l200l iBarth et al.l l2003l~"lDieti-ich et al. 2003 



Freudlin ^etalTl2003l: IJuarez et al.ll2009t) . Although the chemi 
* Based on observations made with the NASA/ESA Hubble Space cal enrichment of the host bkely preda tes the presence of the 
Telescope, obtained at the Space Telescope Science Institute, which is active nucleus dSimon & Hamannl l2010l). winds powered by the 
operated by the Association of Universities for Research in Astronomy, AGN could return this metal-enriched material to the IGM 
Inc., under NASA contract NAS 5-26555. These observations are asso- dFurlanetto & Loebl 12001b ICavaliere etaD l2002t iGermain et"ai1 
ciated with program # 12022. 
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2009; Ho pkins & Elvisll2010t iBarai et alJl201 ll) . Their presence 
and feedback may also have a significant impact on the evolu- 
tion of their host galaxies (Silk & Rees 1998; Scannapieco & Oh 
2004; Gra nato et al1l2004 ; Di Mat teo et alJl2005t iHop kins et al. 
2008; Somerville et al. 2008). While some fraction of these out- 



flows in low-luminosity AGN may not escape their host galaxy , 
at least as me asured in the local universe (iDas et alJ [2005; 
iRuiz et al 2005; Das et al. 2007), the impact of the outflow on 
lower-density portions of the host interstellar medium (ISM) 
could provide a significant transport mechanism for the enrich - 
ment of the surrounding environment (Hopkins & Elvis 2010). 
From constraints provided by the X-ray background, such lower 
luminosity AGN may dominate the population o f active galaxies 
in the early universe dTreister et al.ll2009l 1201 Ol) . 

Nearby AGN provide local analogs that can help us to un- 
derstand the mechanics, energetics, and chemical enrichment 
patterns that may play a significant role in cosmic evolution 
at high redshift. More than half of low-redshift AGN exhibit 
blue-shifted UV or X-ray absorption features indicative of out- 
flowing gas (Crenshaw et alj2003l:lbunn et al.l2007tlCappi et al.1 
2009; Tombesi et al. 2010). Understanding the geometry and 
the location of the outflow relative to the active nucleus is a 
key to making an accurate assessment of the total mass and 
the kinetic luminosity of the outflow. Distance determinations 
are particularly difficult. Using density-sensitive absorption lines 
to establish the gas density, in combination with photoioniza- 
tion models that reproduce the observed relative column den- 
sities can provide precise distance estimates. The se measures 
have r anged from tens o f parsecs in NGC 3 783 dGabel et alJ 
12005b!) and NGC 4151 dKraemer et al.1 120061 if one uses the 
correct critical density fo r Cm*), and up to k i loparsec scales in 
some quasars and AG N dHamannet aUBOOU IScott et alJfeOIMl 
Edmonds "eTaLll201ll). UV and X-ray observations of Mr k 279 
(IScott et al.ll2004t lArav et alj|2007l ICostantini et al. 2007) mea- 
sured absolute abundances in the outflow of a local AGN for 
the first time, and they show evidence for enhanced CNO abun- 
dances that suggest contributions from massive stars and AGB 
stellar winds. 

To improve upon these prior studies, we have conducted 
a multiwavelength campaign of coordinated X-ray, UV, and 
optical observations of th e nearby lum i nous S eyfert 1 galaxy 
Mrk 509 (z=0.034397; iHuchra et all (Il993h). A c omplet e 
overview of the campaign is given by iKaastra et alJ (12011b). 
Mrk 509 is an ideal object for study due t o its high fl ux, mod- 
erate luminosity that rivals that of QSOs (IKopvlov et al. 1974), 
and deep, well structured absorption troughs. The blue-shifted 
absorption indicative of a nuclear outflow has been known since 
the earliest UV spectral observations ( IWu et al. 1980; York et al. 
1984). The outflow is also apparent in the blueshifted [Om] 
emissi on detected acros s the face o f Mrk 509 by Phillips et al. 
(1983). More recently iKriss et all d2000l) and iKraemer etldl 
(2003) observed Mrk 509 at high spectral resolution in the UV 
with the Far Ultraviolet Spectroscopic Explorer (FUSE) and the 
Space Telescope Imaging Spectrograph (STIS) on the Hubble 
Space Telescope (HST), respectively. The STIS observations 
were simultaneous with earlier Chandra X-ray grating obser- 
vations ( Yaqoob et al. 2003) that also detect blue-shifted X-ray 
absorption lines. These prior UV observations show the pres- 
ence of more than seven distinct kinematic components, none of 
which could be unambiguously associated with the absorption 
detected in the X-ray, although both the X-ray and UV absorp- 
tion do overlap in velocity. 

Our newer observations using the Cosmic Origins 
Spectrograph (COS) on HST achieve higher signal-to-noise 



ratios (S/N) and provide insight into long-term changes in the 
absorbing gas. Our COS spectra provide a higher-resolution 
view of the velocity structure in the outflowing gas, and sample 
lower ionization ionic species than are present in the X-ray 
spectra. At high spectral resolution and high S/N, we can use 
velocity-resolved spectral coverage of the N v and C iv doublets 
to solve for the covering fractions and column densities in the 
outflowing gas ( Arav et al. 2007). Using our COS observations 
of Lyo- together with the archival FUSE observations of the 
higher-order Lyman lines, we can provide an absolute reference 
for abundances relative to hydrogen, which cannot be mea- 
sured directly in the X-ray spectra. These higher S/N spectra 
combined with the higher S/N X-ray spectra permit a more 
detailed examination of the links between the UV and the X-ray 
absorbing gas in the outflow. 

This paper describes the COS observations that are part of 
our multiwavelength campaign on Mrk 509, and presents an ini- 
tial interpretation of the results. In §2 we describe our COS ob- 
servations and our data reduction methods, including enhanced 
calibrations for the COS data and corrections for the line-spread 
function. In §3 we describe our analysis of the UV spectra and 
compare them to prior STIS and FUSE observations of Mrk 509. 
Finally, in §4 we discuss the physical implications of our obser- 
vations for the geometry and location of the absorbing gas in 
Mrk 509 and its physical conditions. We end with a summary of 
our conclusions in §5. 



2. Observations and Data Reduction 

As part of our exten sive coordinated campaign on Mrk 509 
(IKaastra et alJ 1201 lbl) . we observed Mrk 509 using the Far- 
Ultraviolet Channel and the medium-resolution gratings of the 
Cosmic Origins Spectrograph (COS) on the Hubble Space 
Telescope (HST). Descriptions of COS and its on-orbit per- 
formance can be found in the COS Instrument Handbook 
(iDixon et al.l 120 10). Our two visits on 2009 December 10 and 
1 1 were simultaneous with the Chandra observations described 
bv lEbrero efail d2011l) . Using gratings G130M and G160M and 
observing through the Primary Science Aperture (PSA), we cov- 
ered the far-ultraviolet spectral range from 1155 A to 1760 A. 
The instrument configurations, times of observation, and inte- 
gration times are summarized in Table Q] For each grating, we 
used only two tilts to avoid placing gaps in spectral ranges of 
interest surrounding the strong emission and absorption features 
in Mrk 509. This resulted in almost complete spectral coverage 
except for a small gap from 1560-1590 A. The two tilts and 
two exposures at each setting resulted in four independent place- 
ments of the grid-wire shadows and other instrumental artifacts 
along the spectrum. This enhanced our ability to flat field these 
high signal-to-noise ratio data. Exposure times were weighted 
more heavily in favor of grating G160M since Mrk 509 is fainter 
across this wavelength range, and the throughput of G160M is 
lower relative to G130M. This gives more uniform signal-to- 
noise ratios across the full spectrum. For G130M we obtained 
a total exposure time of 9470 s; for G160M, the total exposure 
time was 16452 s. 

Our data were processed with the COS calibration pipeline 
v2.1 lb at STScI. This version of the pipeline does not correct for 
flat-field features or the time-dependent sensitivity of COS, and 
there are still some residual anomalies in the wavelength cali- 
bration. To produce the best quality summed data set, we used 
a customized series of procedures to make these corrections as 
described in the following sections. 
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Table 1. COS Observations of Mrk 509 



Data Set Name 


Grating/Tilt 


Date 


Start Time 
(GMT) 


Exposure Time 
(s) 


lbdhOlOlO 


G130M/1309 


2009-12-10 


02:48:40 


1993 


lbdh01020 


G130M/1327 


2009-12-10 


04:07:03 


2742 


lbdh01030 


G160M/1577 


2009-12-10 


05:42:57 


5484 


lbdh01040 


G160M/1589 


2009-12-10 


08:54:42 


2742 


lbdh02010 


G130M/1309 


2009-12-11 


02:46:46 


1993 


lbdh02020 


G130M/1327 


2009-12-11 


04:05:12 


2742 


lbdh02030 


G160M/1577 


2009-12-11 


05:41:05 


2742 


lbdh02040 


G160M/1589 


2009-12-11 


07:16:57 


5484 



2.1. Flat-Field Corrections 

To fiat field the data and correct for instrumental artifacts, we 
used one-dimensional flat-field correcti ons develo ped from COS 
observations of white dwarf standards (Ake et al. 2010). These 
flat fields correspond to each detector segment (A and B) for 
each grating (G130M and G160M) we used. The signal-to-noise 
ratio in these flat fields (in terms of purely Poisson-distributed 
noise) is comparable to that in our data set. To reduce the impact 
of the added noise that would result from dividing by these data 
for our flat-field correction, we first smoothed the flat-fields. It 
may seem counterintuitive to consider smoothing pixel-to-pixel 
flats when using them to correct for pattern noise and artifacts, 
but there are two fundamental reasons why this works. First, 
real artifacts and pattern noise are never sharper than the intrin- 
sic resolution of the detector, which is approximately six pix- 
els. Second, drifts in the grating select mechanism, which are 
corrected during an observation using the TAG-FLASH observ- 
ing mode, effectively smear the pattern noise. Since our obser- 
vations were longer than any of the individual ones used to de- 
velop the one-dimensional flats, the smearing of the pattern noise 
is greater. Smoothing the flats partially corrects for this smear- 
ing. After trying a variety of smoothing widths, we settled on 
smoothing the flats with a Gaussian profile having a 1 -pixel dis- 
persion. To align the flats with the extracted ID spectra from 
each exposure, we cross correlated the appropriate smoothed flat 
with the data, applied an appropriate integer-pixel shift, and then 
divided it out. 

Once the data were flat-field corrected, we aligned the indi- 
vidual exposures by first cross-correlating them with each other 
and applying the closest integer-pixel shift. While this process 
produced a tremendous improvement in the data, some residual 
features remained. We then examined each individual exposure, 
comparing them with each other and with the flat field used to 
correct them. Strong features in the flats that did not divide out 
well were then flagged in the individual exposures so that they 
could be excluded when we did the final combination. Our final 
spectrum for each grating is the exposure-weighted combination 
of the unflagged pixels in each of the individual flat-fielded and 
aligned exposures. Although our data would support a signal-to- 
noise ratio of over 100 per pixel in the Poisson limit, we do not 
do much better than 60: 1 per pixel in our final spectrum. 



2.2. Wavelength Calibration 

While the nominal wavelength calibration of COS over most of 
the spectral range has errors less than 15 kms~', this is not ad- 
equate for our analysis. We wish to align spectral features (both 
intrinsic to Mrk 509 and in intervening material) more accurately 
for comparison in velocity space, and for comparison to prior ob- 
servations with FUSE and STIS. There can also be residual zero- 



point shifts in the wavelength scale due to offsets from the aper- 
ture center which may not have been fully removed in the target 
acquisition process. Therefore, we first adjusted the wavelength- 
scale zero points for our G130M and G160M spectra by cross 
correlating them with the prior STIS observation of Mrk 509. 
For G130M, we used the wavelength region from 1245-1290 
A, which is dominated by the intrinsic absorption features from 
Lyar and N v. For G160M, we used the region from 1590-1615 
A, centered on the intrinsic absorption from C iv. The intrinsic 
absorption features provide numerous sharp edges that result in 
a strong, narrow cross-correlation peak that enhances the preci- 
sion of the alignment. The resulting corrections were < 0. 1 A. 

For the next level of correction, we used the numerous inter- 
stellar absorption lines in our spectra that spanned the full wave- 
length range. As the sight line to Mrk 509 is very complex, each 
interstellar line typically has multiple components . The H I 21- 
cm emission in this direction (Murphy et alJI 19961) has a strong 
flux-weighted peak at +3.1 kms -1 in the local standard of rest 
(LSR), and a weaker, broader peak at — h65 kms -1 (LSR). The 
+3 km s feature is present in all detected interstellar lines, and 
the +65 kms ~' one is present in the stronger lines. In addition, 
there are high-velocity features at -249 and -295 kms visi- 
ble in the strongest lines and in the higher ionization species like 
C iv. We measured the wavelengths of all these interstellar fea- 
tures. The wavelength differences of our measurements relative 
to vacuum wavelengths show a systematic, nearly linear trend 
with wavelength across each detector segment that amount to an 
offset of 0.06 A from one end of each segment to the other. We 
take this as an indication that there is a slight error in the disper- 
sion constant for the COS wavelength calibration. To weight our 
correction by the strongest features, we used only those corre- 
sponding to the +3 and +65 km s components and performed a 
linear fit to the wavelength offsets for each detector segment. We 
then made a linear correction to the wavelength scale, with pivot 
points that have zero correction about 1260 A and 1600 A so that 
we keep our zero points aligned with the STIS wavelength scale. 
After applying this correction, we then find a peak-to-peak vari- 
ation among the velocity offsets of all the interstellar features of 
-10 to +8 kms -1 , with a dispersion of 4.3 kms . 



2.3. Flux Calibration 

The se nsitivity of COS ha s been slowly declining since its instal- 
lation (lOsten et alJ uOlO). We have applied the time-dependen t 
sensitivity corrections to our data as given bv lOsten et alJ (l2010). 
These corrections amount to an adjustment of only 2.5% at the 
short wavelength end of the G130M spectrum, and range up to 
a correction of 5% at the long- wavelength end of G160M. We 
estimate an absolute flux accuracy of 5% at wavelengths above 
1220 A, and 10% at shorter wavelengths dMassa et al.l2010l) . We 
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Fig. 1. Calibrated and merged COS spectrum of Mrk 509. Prominent emission features are labeled above the spectrum. Regions of 
intrinsic absorption are indicated below the spectrum. Geocoronal emission in the center of the galactic Lya absorption trough is 
indicated with an Earth symbol. 



expect that relative fluxes are more precise. Our repeat observa- 
tions in our two visits agree to better than 1%, and compari- 
son between the overlapping wavelength ranges of G130M and 
G160M in our spectra show that they also agree to better than 
1%. 

As an external comparison, and for cross-calibration pur- 
poses with the other elements of our campaign, we compare 
our spectrum to the Swift UVM2 data point from 2010 Dec 10. 
These data are in a broad band centered at 22 31 A, with some 
conta mination from broad Cm] emission (Mehdipou ret all 
2011). To compare our COS spectrum to this point requires 
some ex trapolation. We tak e the 1 140-8090 A spectrum assem- 
bled bv lShangetaH J2005I) from FOS and KPNO 2.1-m spec- 
tra, and normalize this spectrum to the corrected Swift UVM2 
continuum flux at 223 1 A of 5.12 x 10~ 14 ergs cirT 2 s _1 A -1 
dMehdipour et all l201lh . At wavelengths longer than 1220 A, 
our COS spectra agree with the FOS spectrum normalized to the 
corrected Swift UVM2 flux to better than 4% at all wavelengths. 
Below 1210 A, our COS spectrum is brighter by up to 9% above 
the renormalized FOS spectrum. This could be a true difference 
in the spectrum, but it could also be an indication of our residual 
uncertainties in the flux calibration. In the latter case, we conser- 
vatively estimate that our absolute flux below 1220 A is better 
than 10%. 

Figure Q] shows our calibrated and merged COS spectra. We 
have labeled the prominent emission features in the spectrum as 
well as the regions containing the intrinsic absorption features. 
The many galactic interstellar absorption lines in our spectrum 
are not labeled here. In an on-line appendix available only in the 
journal version of this paper, Fig. Al shows full-resolution plots 
of the COS spectrum of Mrk 509. 

2.4. STIS Spectrum 

The STIS spectrum is a 7,600 s observation obtained on 2001- 
04-13 usi ng the echelle E140M grating and the 0.2 x0.2 arc sec 
aperture dKraemer et al.ll2003t lYaqoob et alll2003l Proposal ID 
8877, PI: Yaqoob). The spectrum we use is the one-dimensional 



extracted echelle spectrum direct from the MAST archive with 
all up-to-date calibrations applie d, including correcti ons for scat- 
tered light in the echelle mode ( Valenti et alj|2002h and for the 
echelle blaze corrections ( Alo isi et alj20 07). The STIS spectrum 
of Mrk 509 has a resolution about twice that of the COS spec- 
tra, R=46,000. This spectrum is particularly useful for assessing 
the true resolution of the COS spectrum, and the validity of de- 
convolutions that repair the broad wings of the COS line-spread 
function. 

To illustrate the dramatic improvement in the quality of our 
COS data compared to the prior STIS observation, Figure [2] 
shows a portion of both spectra in the region centered around 
the Lya absorption feature intrinsic to Mrk 509. One sees that 
Mrk 509 was ~60% brighter during our COS observation. In 
addition, the higher throughput of COS also improves our re- 
sulting signal-to-noise ratio by a factor of ~5. However, one can 
see that there are still instrumental features that we must cor- 
rect before doing a detailed analysis of the COS data. Deep, 
saturated absorption features in the STIS spectrum, especially 
the interstellar Si 11 line at 1260.5 A, do not have square, black 
troughs in the COS spectrum. Thi s is due to the br oad wings 
on the line-spread function (LSF) (iGhavamian e t al. 2009) that 
redistribute light from the bright continuum into cores of the ab- 
sorption lines. To extract the maximum information from our 
spectrum, we must correct for these broad wings on the LSF. 



2.5. Deconvolving the COS Spectrum 

The broad wings on the COS LSF are caused by mid-frequency 
polishing errors on the HST primary & secondary mirrors 
(Ghavamian et al. 2009). Additional scattering by the HST mir- 
ror system and internal to C OS makes these wings more exten- 
sive than originally thought dKrissll201lT> . These errors must be 
removed in order to obtain accurate measurements of the depths 
of the absorption lines in the COS spectrum. On scales less than 
about 50 kms~ , significant light can leak into the absorption 
line troughs from adjacent continuum regions. As a result, in- 
terstellar absorption lines expected to be saturated are not black. 
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Fig. 2. COS spectrum of Mrk 509 in the Lya region (black) com- 
pared to the STIS spectrum (red). 



Fig. 3. Comparison of the deconvolved COS spectrum in the 
spectral region surrounding the Lya absorption trough (black) to 
the original COS spectrum (green) and the STIS spectrum (red). 



This is quite obvious when one compares the interstellar lines in 
the COS spectrum in Fig. [2] to the prior spectrum obtained with 
STIS. 

The high signal-to-noise ratio of our data permits us to 
do a fairly good job of deconvolving the line-spread function 
from our data. We use the Lucy-Richardson algorithm to decon- 
volve the data as implemented in the STSDAS routine "lucy" 
in the analysis.restore package. We deconvolve the spectrum in 
50 A intervals using the updated versions of the wavelength- 
dependent LSFs from the COS web site at STScI dKrissll201ll) . 
The deconvolution converges well after ~20 iterations. A good 
illustration of the success of our deconvolution is provided by 
direct comparisons to the higher-resolution STIS spectrum of 
Mrk 509. Figure [3] shows the normalized version of the origi- 
nal COS spectrum in the Lya region, the deconvolved spectrum, 
and the STIS spectrum, all plotted together. The deconvolution 
deepens the narrow interstellar absorption features so that they 
match the higher-resolution STIS spectra well. Broader, satu- 
rated interstellar lines like Si n A 1260.5 have square, black bot- 
toms. Overall, the deconvolved, high S/N COS spectrum looks 
like an excellent fit to the STIS spectrum. One disadvantage to 
the deconvolution, however, is that it amplifies noise in the con- 
tinuum. This can obscure weak absorption features. Therefore, 
in our analysis we will adopt the approach of using the origi- 
nal spectrum to identify significant features, but use the decon- 
volved spectrum to make accurate measurements of the strength 
of known features so that the effects of the LSF are fully cor- 
rected. 



2.6. FUSE Spectra 

FUSE (lMoosetal.ll2000t) observed Mrk 509 on two different 
occasions during the mission. The first FUSE spectrum, ob- 
tained in two closely spaced observations on 1999-1 1-02 (19,355 
s) and 1999-1 1-06 ( 32,457 s), were analyzed and discussed 
bv lKrissetaT](l2000l) . The second FUSE visit (Observation ID 
P1080601) was a 62,100 s observation obtained on 2000-09-05. 
We retrieved these data from the MAST archive and reprocessed 
them with a customized version of the FUSE pipeline, using only 
the night portions of the observation in order to minimize air- 



glow and background. Our version of the processing also low- 
ers the background noise by screening out more of the lowest 
pulse height channels. We cross correlate the overlapping wave- 
length sections of each channel to align the wavelength scales, 
and renormalize all the fluxes to match channel LiFIA, which 
was aligned with the fine guidance sensor for both observations. 
We correct for the "worm" in segment LiF2A by taking a ra- 
tio to the matching portion of the spectrum in LiFIB, fitting a 
spline, and applying that correction to LiF2A. We then combine 
all channels in overlapping wavelength regions weighted by the 
S/N of the individual pixels. Although Mrk 509 was fainter in 
2000 compared to the first observation, channel alignment was 
better, and the overall S/N is higher. Mrk 509 was fainter in both 
FUSE observations than in our current COS observation by 34% 
to 55%. 



3. Data Analysis 

3. 1 . Continuum and Emission Lines 

Our first step in analyzing the absorption features in Mrk 509 
is to fit the line and continuum emission. We use the 
specfit task in the contrib package of IRAF for our fits 
(iKrissI Il994l) . For the continuum, we use a reddened power 
law (F A = Fiooo (A/1000A) _a ) with extinction of E(B- 
V) = 0.057 (as ob tained from NED via the prescription of 
Schlegel et al. 1998& a nd var iation with wavelength as speci- 
fied by Card elli et al.l d!989l) with a ratio of selective to to- 
tal extinction of Ry = 3.1. Our best-fit power law has F^ = 
2.48xl0~ 13 (A/l000k)- lJ4 ergs cnT 2 s" 1 A" 1 . A summary of the 
power law parameters and observed continuum fluxes at 1 175 A 
for the COS observation and the prior STIS and FUSE observa- 
tions is given in Table [2] 

For the strongest emission lines, our model requires four 
Gaussian components to give a good fit. For O vi, Lya, and 
Civ, we use a broad base, two moderately broad components 
that comprise most of the line profile, and a narrow-line compo- 
nent of full-width half-maximum (FWHM) of 300 kms ~ fixed 
to match t he width of the nar row [O in] emission in the visible 
spectrum (Phillips et al.lll983D . For the O vi and Civ doublets, 
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Table 2. Continuum Fits to Observations of Mrk 509 



Observatory 


Date 




Fiooo 




a 




F(1175 A) 






(io- 


13 -2 - 

ergs cm s 


'A- 1 ) 




(io- 


13 ergscirT 2 s _1 A~') 


FUSE 


1999-11-06 




1.58 




0.41 




0.607 


FUSE 


2000-09-05 




0.88 




1.16 




0.416 


HST/STIS 


2001-04-13 




1.44 




1.74 




0.538 


HST/COS 


2009-12-10 




2.48 




1.74 




0.918 



we include a component for each line of the doublet with the 
relative wavelengths fixed at the ratio of the laboratory values, 
and an optically thin, 2:1 ratio for the blue to the red flux com- 
ponent of the doublet. Weaker lines in the spectrum (e.g., C in, 
N v, Si iv) require only a moderately broad component. The only 
other visible narrow line in the far-UV spectral range is a He n 
/11640 line. The panels in Fig. [4] show the full emission model 
overplotted on the COS spectrum of Mrk 509 and the individual 
emission components for the regions surrounding the Lya, N v, 
and Civ emission lines. Tables [3] H E and [6] give the best fit- 
ting parameters for the fitted emission lines in the COS, STIS, 
and FUSE spectra. We note that the narrow-line components are 
weak, and they are poorly constrained. They do improve the fit, 
but since they have such a small contribution to the unabsorbed 
spectrum, they ultimately have little impact on the properties de- 
rived for the absorption lines. 

3.2. Foreground absorption features 

Given our description of the emission lines and continuum in 
our spectrum of Mrk 509, we now identify and fit all signifi- 
cant absorption features in the observed spectrum. Most of these 
features are foreground absorption due to our own interstellar 
medium, and they will be discussed in detail by Pinto (in preps.). 
We fit the galactic and intergalactic absorption features using a 
Gaussian profile in optical depth. Many of the strongest inter- 
stellar lines are blends of features at several different velocities. 
Up to seven individual components are detected in some ISM 
lines. In cases where we cannot definitively separate these com- 
ponents, we fix their relative velocities at the mean of the ve- 
locities observed in cleanly resolved components of other inter- 
stellar lines. Our measured velocities and equivalent widths for 
the identified features are given in Table [7] The tabulated error 
bars are purely statistical. For the velocities, at a resolution of 15 
km s _1 and S/N > 50, the statistical errors are expected to be < 1 
kms _I . In such cases, we have rounded the values up to 1. 

We also find three foreground absorption features that we 
identify as Lye absorption in the intergalactic medium as shown 
in Table [8] One of these features was previously identified by 
iPenton et al.l d2000). The second feature at 1221.08 A is simi- 
lar to this, but much weaker. The third feature at 1239.03 A is 
badly blended with foreground galactic N v /11238. Given the 
observed equivalent width of the red component of N v A 1242, 
this blend is too strong to be simply foreground N v absorption. 
We attribute the excess (and the slight velocity offset) to an in- 
tergalactic Lyff absorber. We observe no other absorption lines 
in our spectrum associated with these foreground intergalactic 
systems. 



3.3. The Intrinsic Absorption Features 

The COS spectrum of Mrk 509 shows a wealth of detail in the 
intrinsic absorption features. The sev en major components origi- 
nally identified by Kriss et al.l d2000t) in the FUSE spectrum were 
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Fig. 4. Fit to the broad lines and continuum of Mrk 509 (red line) 
is plotted on the observed spectrum (black). Green lines show the 
individual emission-line components in the fit. The blue line is 
the continuum. The first panel shows the Lya and N v region. 
The second panel shows the C iv region. 



expanded to eight bv lKraemer et al.l (120031) based on their STIS 
spectrum. Our COS spectrum shows additional subtle inflections 
indicating that the absorption troughs are a quite complex blend 
of perhaps even more components than we can readily resolve. 
Figure [5] shows the intrinsic absorption features detected in our 
COS spectrum as a function of velocity with the components 
labeled. For reference, we also include the O vi, Ly/3 and Lyy 
regions of the spectrum from the FUSE observation of Mrk 509 
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Table 7. Interstellar Absorption Features 
in the COS Spectrum of Mrk 509 



Table 7. Interstellar Absorption Features 
in the COS Spectrum of Mrk 509 (cont.) 



Feature 


^0 


EW 


Vlsr 




(A) 


(mA) 


(km s" 1 ) 


Ci 


1157.91 


15 ±3 


-9+1 


S in 


1190.21 


150 ±18 


+2+1 


Sin 


1190.42 


75 ±7 


-48 + 2 


Sin 


1190.42 


306 ± 13 


+7+1 


Sin 


1190.42 


360 ± 34 


+65+ 1 


Sin 


1190.42 


11 ±2 


+143 ± 1 


Ci 


1192.22 


10 ±2 


-6+1 


Sin 


1193.29 


56 ±4 


-52+ 1 


Sin 


1193.29 


329 ± 18 


+ 11 + 1 


Sin 


1193.29 


345 ± 28 


+72+ 1 


Sin 


1193.29 


12±2 


+ 132+1 


Mnn 


1197.18 


24 ±3 


+5 + 1 


Mnn 


1199.39 


45 ±6 


+3+1 


Mnn 


1201.12 


14 ±2 


+4+1 


Ni 


1199.55 


193 ±5 


+9+1 


Ni 


1199.55 


171 ±4 


+74+1 


Ni 


1200.22 


173 ±5 


+8+1 


Ni 


1200.22 


149 + 3 


+70+1 


Ni 


1200.71 


175 ±5 


+6+1 


Ni 


1200.71 


115±3 


+69+1 


Sim 


1206.50 


202 ±4 


-289 ± 1 


Sim 


1206.50 


59 ±3 


-227 ± 1 


Sim 


1206.50 


11 ± 1 


-63+ 1 


Sim 


1206.50 


417 + 37 


+6+1 


Sim 


1206.50 


477 ± 84 


+66+ 1 


Sim 


1206.50 


26 + 2 


+140+1 


Nv 


1238.82 


11+2 


-292 ± 1 


Nv 


1238.82 


7 + 2 


-236 ± 1 


Nv 


1238.82 


37 + 2 


+ 11 + 1 


Nv 


1238.82 


10 + 2 


+69+ 1 


Nv 


1242.80 


6 + 2 


-294 ± 1 


Nv 


1242.80 


4 + 2 


-238 ± 1 


Nv 


1242.80 


19 + 2 


+ 10+1 


Nv 


1242.80 


5 + 2 


+68+ 1 


Mgn 


1239.93 


27 + 2 


+ 10+1 


Mgn 


1240.39 


16 + 2 


+8+1 


Sn 


1250.58 


117 + 3 


+4+1 


Sn 


1250.58 


34 + 2 


+60+ 1 


Sn 


1253.81 


138 + 3 


+8+1 


Sn 


1253.81 


72+1 


+65 + 1 


Sn 


1259.52 


151+2 


+4+1 


Sn 


1259.52 


85 + 1 


+64+1 


Sn 


1259.52 


153+ 1 


+4+1 


Sn 


1259.52 


86+14 


+64+1 


Sin 


1260.42 


14 + 2 


-291 ± 1 


Sin 


1260.42 


6 + 3 


-96+ 1 


Sin 


1260.42 


51+6 


-61 + 1 


Sin 


1260.42 


442 ± 15 


+ 11 + 1 


Sin 


1260.42 


341 + 11 


+75+ 1 


Sin 


1260.42 


31+4 


+ 134+1 


Ci 


1277.25 


51 + 32 


+5 + 1 


Pii 


1301.87 


41+3 


+3+1 


Oi 


1302.17 


16 + 6 


-292 ± 1 


Oi 


1302.17 


6 + 6 


-97+ 1 


Oi 


1302.17 


2 + 6 


-64+ 1 


Oi 


1302.17 


307 + 9 


+7+1 


Oi 


1302.17 


277 ± 21 


+74+1 


Oi 


1302.17 


6 + 5 


+ 131 + 1 


Sin 


1304.37 


16 + 6 


-290 ± 1 


Sin 


1304.37 


6 + 6 


-95+ 1 


Sin 


1304.37 


8 + 6 


-39+ 13 


Sin 


1304.37 


284 ± 1 1 


+9+1 


Sin 


1304.37 


251 ± 10 


+66+ 1 


Sin 


1304.37 


15 + 6 


+ 133+ 1 


Ci 


1328.83 


27+4 


-1 + 2 


Cii 


1334.53 


76 + 4 


-293 ± 1 


Cn 


1334.53 


23+4 


-223 ± 1 



Cn 


1334.53 


65 + 4 


-59 + 1 


Cii 


1334.53 


496 ± 37 


+ 11 + 1 


Cn* 


1335.50 


444 ± 21 


-144+1 


Cn* 


1335.50 


150 + 4 


+49+1 


Cn* 


1335.50 


105+4 


+ 105+1 


Nin 


1317.22 


33 + 8 


+3 + 3 


Nin 


1317.22 


16 + 3 


+64+1 


Nin 


1370.13 


41 + 3 


+4+ 1 


Si iv 


1393.76 


95 + 4 


-298 ± 1 


Si iv 


1393.76 


16 + 3 


-246 ± 1 


Si iv 


1393.76 


19 + 3 


-63 + 1 


Si iv 


1393.76 


227 ± 14 


+8 + 2 


Si iv 


1393.76 


165+ 18 


+64 + 3 


Si iv 


1393.76 


7 + 5 


+ 131 + 1 


Si iv 


1402.77 


53 + 3 


-297 ± 1 


Si iv 


1402.77 


11 + 3 


-275 ± 3 


Si iv 


1402.77 


7 + 3 


-64+1 


Si iv 


1402.77 


144 + 7 


+9+ 1 


Si iv 


1402.77 


103 + 8 


+65 + 2 


Si iv 


1402.77 


2+1 


+ 131 + 1 


Nin 


1454.84 


48 + 5 


-6+ 1 


Siiv a 


1393.76 


89 + 4 


-298 ± 1 


Si iv 


1393.76 


14 + 3 


-252 ± 3 


Si iv 


1393.76 


19 + 3 


-67+1 


Si iv 


1393.76 


223 + 8 


+4+ 1 


Si iv 


1393.76 


160 + 8 


+60+1 


Si iv 


1393.76 


20 + 4 


+ 126+1 


Si iv 


1402.77 


65 + 2 


-297 ± 1 


Si iv 


1402.77 


14 + 2 


-250 ± 3 


Si iv 


1402.77 


7 + 2 


-62+ 1 


Si iv 


1402.77 


157 + 7 


+9+ 1 


Si iv 


1402.77 


116 + 8 


+69+1 


Si iv 


1402.77 


2 + 2 


+ 131 + 1 


Nin 


1454.84 


24 + 2 


-2 + 2 


Nin 


1454.84 


11 + 2 


+56+ 1 


Sin 


1526.71 


25 + 2 


-73+ 1 


Sin 


1526.71 


408 ± 1 1 


+0+ 1 


Sin 


1526.71 


340+ 11 


+59+1 


Civ 


1548.19 


245 + 5 


-298 ± 1 


Civ 


1548.19 


78 + 5 


-245 ± 1 


Civ 


1548.19 


19+3 


-125 + 3 


Civ 


1548.19 


22 + 4 


-63+ 1 


Civ 


1548.19 


336 + 8 


+8+ 1 


Civ 


1548.19 


293 + 9 


+66+1 


Civ 


1548.19 


18 + 4 


+ 130+1 


Civ 


1550.77 


170 + 4 


-298 ± 1 


Civ 


1550.77 


52 + 4 


-248 ± 1 


Civ 


1550.77 


1 + 14 


-269 + 11 


Civ 


1550.77 


10 + 4 


-64+1 


Civ 


1550.77 


251+7 


+5+ 1 


Civ 


1550.77 


178 + 7 


+65 + 1 


Civ 


1550.77 


11+4 


+ 129+1 


Ci 


1560.31 


45 + 4 


-2+ 1 


Fen 


1608.45 


22 + 2 


-73 + 1 


Fen 


1608.45 


217 + 5 


+2+ 1 


Fen 


1608.45 


204 + 5 


+59+1 


Fen 


1608.45 


19 + 3 


+91 + 1 


Fen 


1611.20 


36 + 4 


+2+2 


Ci 


1656.93 


72 + 5 


-3+ 1 


Ci 


1656.93 


15 + 3 


+59+1 


Alii 


1670.79 


33 + 3 


-64 + 8 


Alii 


1670.79 


420 ± 30 


+3+ 1 


Alii 


1670.79 


379 ± 28 


+56 + 2 


Alii 


1670.79 


23 + 5 


+ 121+4 


Nin 


1709.60 


63 + 12 


+ 11+4 


Nin 


1709.60 


28 + 6 


+69 + 4 


Nin 


1741.55 


50+10 


+8 + 6 


Nin 


1741.55 


20 + 9 


+66 + 6 


Nin 


1751.91 


36+ 13 


+3 + 5 


Nin 


1751.91 


17+10 


+61+5 



"The prior entries for Si iv /i/il393,1402 and 
Ni ii /U454 were for features in the G130M 
spectrum. These subsequent entries are for 
features detected in the G160M spectrum. 



G. A. Kriss et al.: Multiwavelength campaign on Mrk 509 



Table 8. Intergalactic Absorption Features 
in the COS Spectrum of Mrk 509 



-400 -200 



200 



Feature 



(A) 



EW 

(mA) 



(kms- 1 ) 



Hi 
Hi 
Hi 



1215.67 
1215.67 
1215.67 



31 i 

200: 

72i 



4-1345 ±2 
+2561 ± 1 

4-5772 ± 2 



obtained in 2000. The numbering scheme follows that originated 
by iRriss et al. (2000), but includes the additional component 4' 
added bv lKraemer et al.l (120031) . However, instead of continuing 
to designate additional components with a "prime", we use des- 
ignations of "a" or "b" for the additional features, with the con- 
vention that an "a" subcomponent is blueshifted from the main 
feature, and a " b" subcomponent is redshifted. Component 4' in 
iKraemer et al.l (J2003) becomes 4a here. 

To fit these features with a simple Gaussian decomposi- 
tion of the profiles as done by Kriss et alj (2000) clearly has 
limitations built in by the choice of the number and place- 
ment of components. Given the complexity of the absorp- 
tion troughs, a more model-independent approach is preferred. 
Several previous studies of AGN outflows have used absorption- 
line doublets to measure the optical dept h and covering frac- 
tion in discrete, velocity resolved bins (e.g ., Hamann etal.ll997 



1999; 



2003 



iBarlow & Sareenllll997): lAray et alJ ll999J: iGangulv et all 
NGC 55 48: lAravetalJ 120021: NGC 3783: iGabelet all . _ 
Mrk 279: IScott et alJl2004l iGabel et aDl2005al) . We will apply 
this method to our spectrum of Mrk 509 in a subsequent pub- 
lication (Arav, in preps.). In this paper, however, we obtain a 
preliminary assessment of the column densities and covering 
fractions of the intrinsic absorption components by fitting the 
features with Gaussian profiles in optical depth. This enables a 
more straightforward comparison of our current results to the 
prior FUSE and STIS observations. Our model also allows for 
partial covering by each component, as well as a scattered light 
component that is a wavelength-independent fraction of the un- 
derlying emission components. We note that in the case of a sin- 
gle absorbing region with a single partial covering fraction (de- 
noted f c ) that this is observationally indistinguishable from hav- 
ing a scattered light component at a level of ( 1 - f c ). However, in 
the case where there is more than one absorbing component that 
overlaps in velocity with another absorbing component, there is 
a significant difference. Assuming that the absorbers are statis- 
tically independent, in the region of velocity /wavelength space 
that they overlap, the transmission of light from the source will 
be less than for either component individually, being the prod- 
uct (1 - / c i)(l - fa)- This produces extra dips in the modeled 
spectrum in the regions of overlap as illustrated in Fig. [6] This is 
not a qualitatively accurate representation of the profiles of the 
absorption troughs in our spectrum, which tend to show more 
smoothly varying profiles. Characterizing the residual light at 
the bottoms of our absorption troughs with a scattered compo- 
nent that is independent of wavelength provides a better match 
to the absorption profiles that we observe. For the N v, Si iv, 
and C iv doublets, we link the relative wavelengths at the ratio of 
their vacuum wavelengths, and we link the relative optical depths 
of the red and blue components at a value of 0.5. Since the Lya 
absorption appears to be heavily saturated, as noted previously 
by K raemer et al.l (120031) . we fix the covering fraction for each 
of the Lya components at the average of the value determined 
for the corresponding components in the C iv and N v profiles. 




-400 -200 



200 



Velocity (km s -1 ) 

Fig. 5. Intrinsic absorption features in the COS and FUSE spec- 
tra of Mrk 509. Normalized relative fluxes are plotted as a func- 
tion of velocity relative to the systemic redshift of z — 0.034397. 
The top panel shows the O vi doublet from the FUSE observa- 
tion in 2000 September with the red side of the doublet as a red 
line and the blue side of the doublet as a blue line. The mid- 
dle four panels show the N v, C iv, Si iv and Lya absorption 
troughs from our COS observation. The bottom panel shows the 
Lyy (red) and Ly/? (black) absorption troughs from the FUSE 
observation. Dotted vertical lines indicate the velocities of iden- 
tified absorption components. The thin black vertical line shows 
v = 0. 

Our best-fit resu l ts are given in Table [9] Note that for Si iv, 
as in IKraemer et al.l (120031) . we detect only component 4. Also, 
as in Kraemer et al. (2003), no lower ionization lines (e.g., Si in 
/U206, Si n ,11260, or Cn ,11335) are detected. As we noted 
in §3.1, the narrow emission-line components in our emission 
model have little impact on the measured properties of the ab- 
sorption lines. Even if we eliminate the narrow-line emission 
from our fits, the measured columns densities change by < 1%. 

Using the additional components that we have identified in 
our COS spectra, we have re-analyzed the archival FUSE and 
STIS spectra of Mrk 509 in order to make a direct comparison 
among the spectra. While iKriss et alj (J2000) have done this for 
the first FUSE observation, similar measurements for the sec- 
ond observation are useful since they illustrate how the absorbers 
may change in response to changes in the continuum strength. 

Our results from fitting the STIS spectrum are given in 
Table [10] The column densities for the individual components 
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1257 1257.5 1258 1258.5 

Wavelength (A) 

Fig. 6. Fit to a portion of the Lya absorption trough in the COS 
spectrum of Mrk 509. The data are shown by the solid black his- 
togram. The red dashed line shows the best fit model using the 
scattered-light/single partial covering fraction described above 
for f c - 0.064. The black dashed line shows the emission model 
for this region. The absorption in this model due to components 
4a (shown in green), 4 (shown in blue) and 5 (shown in green) 
are plotted individually. The solid red line shows the resulting 
model if these absorbing components each have independent 
partial covering factors f c ^ a - f c ,4 - fc,5 - 0.064. The dips 
in the region of overlap between each component are prominent 
at wavelengths of ~ 1257.3 A and 1257.6 A. 



agree well (and within the errors) with those measured by 
Kraemeret al. (2003) except for the blended components com- 
prising features 2 and 3, and 6 and 7. The apportionment among 
the various subcomponents accounts for these differences in our 
fit. We also find a higher column density in the strongest fea- 
ture, number 4, largely because o ur width for that com ponent is 
narrower than that measured by Rraemer et al. (2003). 

Results from fits to Ly/3, Lyy and O vi in the 1999 FUSE 
spectrum are in Table QT| and results for the 2000 FUSE ob- 
servation are in Table [12] Since the Lya absorption trough in 
Mrk 509 is heavily saturated, the higher-order Lyman lines with 
good S/N in the FUSE spectrum are particularly useful in de- 
riving reliable column densities and covering factors for H i. 
However, given potential variability, these values must be used 
with caution when combined with the data from our COS obser- 
vations. 

In carrying out our fits we found that the well defined fea- 
tures in N v and C iv in the COS spectrum did not always match 
well with the absorption profiles in Lya and in O vi. One can 
see this in Fig. [5] where the Lya and O vi profiles vary more 
smoothly. For our best-fit results in Tables |9T[T2"1 one can see 
that some components, such as la, are significant only in Lya-, 
while others, such as lb, are strongest in O vi. There are also 
small velocity shifts for other components when comparing N v 
and C iv to Lya and O vi. This illustrates the limitations of our 
Gaussian decomposition, and in §4.2 of the Discussion we ad- 
dress the probable presence of multiple phases of gas at each 
velocity in the intrinsic absorption troughs. 



3.4. Comparison among the COS, STIS and FUSE Spectra 

As one can see from the fits to the continua of each of the 
observations, we have sampled the emission from Mrk 509 
over a range of about a factor of two in overall luminos- 
ity. The FUSE observation in September 2000 represents the 
faintest state (0.42 x 10~ 13 ergscrrr 2 s~' A -1 at 1175 A), and 
the COS observation reported here is the brightest (0.92 x 
10~ 13 ergscm^s -1 A -1 ). 

Our COS spectrum of Mrk 509 captures the active nucleus 
in a significantly brighter state compared to STIS. Overall, the 
continuum is ~80% brighter. Despite the difference in bright- 
ness, the continuum shape does not noticeably change between 
the STIS and the COS observation. The broad emission lines 
are also brighter, but not by as much — only 37%. Despite the 
difference in the continuum level, there are only subtle differ- 
ences in the absorption-line troughs. The much higher S/N in 
our COS spectrum reveals some interesting differences from the 
original STIS spectrum. We show this by comparing the de- 
convolved COS absorption-line profiles to the STIS profiles in 
Fig- [7] The biggest differences are in the Lya absorption profile. 
Significantly less absorption is present in the most-blue-shifted 
portions of the trough (essentially components 1— lb). The most 
striking difference, however, is the extra light that seems to fill 
the bottom of the Lya absorption troughs in the COS spectrum, 
even in the deconvolved version. Comparable differences in the 
N v and C iv troughs are not readily visible in Fig. [7j so we have 
carefully examined the COS data to rule out any instrumental 
effects that might be causing this. 

First, we can rule out instrumental scattered light, or resid- 
ual effects of the line spread function. Since different grating tilts 
illuminate the detector slightly differently, we compared the sep- 
arate spectra from each grating tilt and each visit. The extra light 
in the Lya trough is visible at both grating tilts, 1309 and 1327, 
in both visits of our observations. One can readily see in the de- 
convolved COS spectrum in Fig.[3]that the saturated galactic Si n 
A 1260 absorption line is fully corrected to a black trough that is 
a good match to the STIS spectrum. The much deeper and wider 
absorption trough at galactic Lya at 1216 A is also black, even 
in the spectrum that has not been deconvolved if one looks at 
regions that are more than a resolution element removed from 
the geocoronal Lya line at the center of the trough. (The peak 
of the geocoronal Lya emission is far brighter than the peak of 
Lya in Mrk 509, yet the flux drops quickly by over two orders of 
magnitude from the peak down to the trough of the galactic Lya 
absorption.) 

Another possibility is that extended emission in Mrk 509, 
in particular, extended narrow Lya emission might be filling in 
the absorption trough. The STIS observation used the 0.2x0.2 
arc sec aperture, while the COS entrance aperture is 2.5" in di- 
ameter. Extended emission that is excluded by the STIS aper- 
ture could be entering the COS aperture. We examined the 2D 
spectra, both visually and as summed projections showing the 
spatial profile, and there is no difference in the spatial extent of 
the light in the Lya troughs compared to the continuum. The 
cross-dispersion profile of the continuum spectrum also shows 
the nominal width and shape of a point source in the COS aper- 
ture. In contrast, the spatial extent of the geocoronal Lya emis- 
sion, which fills the aperture, is quite obvious. Of course it is 
quite bright, but even scaled down it is easily seen to be ex- 
tended. We conclude that there is no evidence for extended line 
emission visible in the COS aperture. 

Looking closely at the STIS spectrum, one sees that the Lya 
absorption is not black at all portions of the trough. In fact, 
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Fig. 7. Comparison of spectral features in the COS (black) and 
STIS (red) spectra of Mrk 509. Normalized relative fluxes are 
plotted as a function of velocity relative to the systemic redshift 
of z - 0.034397. The individual panels (from top to bottom) 
show Lya, N v ,11238, N v ,11242, C iv ,11548, and C iv ,11550. 



iKraemer et al.l (120031) required components with partial cover- 
ing to fit the intrinsic absorption lines in their spectrum (C iv and 
N v as well as Lya), as did Kriss et al. (2000) in their FUSE 
spectrum. If one compares the partial covering fractions among 
all three data sets (FUSE, STIS, and COS), one sees that COS 
and FUSE require similar covering fractions, and STIS some- 
what less. The to tal depth of the O vi absorption troughs in 
iKriss et al . (2000) are almost directly comparable to the most 
highly saturated portions of the Lyor trough in the COS spec- 
trum, with the deepest portions having ~ 92% coverage (includ- 
ing both the scattered light components and the partial covering 
factors). Note that both the COS and FUSE spectra were ob- 
tained through large apertures (2.5" diameter and 30" x 30", re- 
spectively), suggesting that an extended emission component (or 
light from an extended scattering region) is filling in the absorp- 
tion troughs. From our inspection of the COS two-dimensional 
images described above, since we do not resolve this component 
in the cross-dispersion direction, we conclude that such a region 
must be larger than the STIS aperture size, but smaller than the 
nominal COS cross-dispersion spatial resolution, which is 0.5" 
at its best at a wavelength of 1600 A (Ghavamian et al. 2010). 
Comparing the best-fit scattered light components from our fits 
to the COS and STIS spectra, 0.05 vs. 0.01, if we assume that 



Velocity (km s -1 ) 

Fig. 8. Comparison of spectral features in the FUSE spectra of 
Mrk 509 from 1999 November (blue) and 2000 September (red). 
Normalized relative fluxes are plotted as a function of velocity 
relative to the systemic redshift of z — 0.034397. The individual 
panels (from top to bottom) show O vi ,11032, O vi ,11037, Ly/?, 
and Lyy. Foreground Galactic ISM features are marked. Error 
vectors for each spectrum are shown as thin lines in the same 
color as the data. 



this light is coming from a region of uniform surface brightness, 
the total extent of this region would be V5 times the size of the 
STIS aperture, or roughly 0.45" in diameter, which is consistent 
with the limits in size we set from the lack of any detectable 
spatial extent in the COS spectrum. 

This scattering region would then have a radius of ~ 150 
pc. If it occupied a biconical region with a covering fraction of 
AQ/4tt = 0.25, following Kroli k & Begelmanl dl986h . the re- 
quired Thomson optical depth T e to produce the scattered flux 
would be T e = 0.05/(AQ/4tt) = 0.20, where 0.05 is the fraction 
of uniformly scattered light in the troughs of the COS and FUSE 
spectra. This gives a total column density of 3.0 x 10 23 cm -2 for 
the hot gas. If it uniformly fills this volume, it has a density of 
1900 cirT 3 and a total mass of 2.0 x 10 8 M @ . While this is larger 
than the reflecting region in NGC 1068, its properties are not 
unreasonable for a volume filled by the outflow from the nuclear 
region of Mrk 509. 

In addition to this scattered light component, we see dif- 
ferences in the depths of several absorption components in the 
Lyor and N v profiles. Examining Fig. [7] once again, one can 
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see that components 1 and lb in the COS Lya spectrum are not 
as deep as in the STIS spectrum. In contrast, component 1 in 
the N v absorption profile is noticeably shallower in the STIS 
spectrum compared to the COS spectrum. More quantitatively, 
we can compare the column densities in Tables l9l and [TUl From 
the FUSE Lyy and Ly/3 profiles in the velocity range of com- 
ponents lb, 1, and la we can see that this region of velocity 
space is optically thin, even in the Lya profile. Therefore, the 
change in the summed column density from these three compo- 
nents from 88x1 12 cirT 2 to57xl0 12 cirT 2 is a decrease of 35%. 
For N v, the column density increased from 1 12 x 10 12 crrT 2 to 
135 x 10 12 crrT 2 . If these changes are due to a change in the ion- 
ization state of the absorbing gas, it indicates that the gas must 
have been in a lower ionization state in 2000, at an ionization 
parameter somewhat below the peak for the ionization state of 
Nv. 

Changes in absorption are also apparent when one compares 
the two FUSE observations, obtained roughly 10 months apart. 
In Fig. [8] we compare the absorption profiles in velocity space 
for the O vi doublet and for Ly/3 and Lyy. (One must be care- 
ful in comparing these spectra since portions of each line profile 
are contaminated by foreground Galactic absorption.) We have 
marked the contaminating features in each panel of Fig. [8] O vi 
is very heavily saturated at almost all velocities; we see no sig- 
nificant differences visually in Fig. [8] The fits in Tables QT| and 
IT21 also show similar column densities, but, given the level of 
saturation in O vi, these are not reliable comparisons. In com- 
paring the Ly/3 and Lyy profiles, however, one can see signifi- 
cant changes in the red side of the absorption troughs, especially 
for components 5, 6a, 6, 6b. There are also slight differences in 
the region of component 1, as we saw in comparing the COS 
and STIS observations. Doing a quantitative bin-by-bin compar- 
ison in velocity space over the range of components la, 1, and 
lb, and 6a, 6, and 6b, we find that changes in the component 1 
region, which primarily show up in Lyy are significant at only 
the 1-cr confidence level. For components 5, 6a, 6, and 6b, how- 
ever, the changes in Ly/3 and Lyy are significant at > 3cr confi- 
dence. Summing over all of these components for H i in Tables 
fm and [12] we see that the total column of H i increased from 
1.16 x 10 15 cnr 2 to 1.30 x 10 15 cnr 2 from 1999 November to 
2000 September. If due to a change in ionization, this is consis- 
tent with the observed decrease in continuum flux between the 
two observations. 



4. Discussion 

4. 1. The Environment of the Active Nucleus in Mrk 509 

The velocities for the absorbing gas in Mrk 509 span a range 
from negative to positive that is atypical for AGN. When UV 
absorption lines are detected in low-redshift AG N, they ar e 
generally blue shifted (Crenshaw et al. 2003; Dunn et al. 2007). 
This is also true of higher r edshift quasar s dWeymann et al.l 
1979t iGanguly et alj |200U iMisawa et al.l |2007|), although 
Wevmann et alJ (ll979l) did suggest that systems with red-shifted 



absorption lines could be associated with neighboring galax- 
ies in clusters. However, given the significant blue shifts of 
the UV emission lines often used to define the quasar redshift, 
such redshifted systems could well be blueshifte d relative to the 
systemic velocities of the host galaxies. In fact, Gangul y et al.l 
finds that the distribution of associated narrow-line ab- 




Fig.9. HST WFPC2 image of Mrk 509 through filter F547M. 
The directions North and East are marked with arrows of length 
2.5". The black circle shows the placement of the 2.5" diam- 
eter COS aperture used for our observations. Lines indicating 
the position angles of various features associated with the struc- 
ture of Mrk 509 are indicated. The small knot at the end of the 
linear structure visible in the F547M image directly above the 
nucleus lies at a position angle (PA, measured East of North) 
of -42° from the nucleus. The minor axis (PA= 1 35°) of the 
low-ionization gas disk is from Phillips et al.l (Il983l) . The PA of 
the polarized light in the core of the Ha emission line (60°), the 
broa d wings of Ha (168°) and the optical continuum (144°) are 
from You ng etalJ dl999|). The radio structu res (PA= -65° and 
165°) are from Singh & Westergaardl d!992l) . Note that extended 
radio emission to the NW rotates from PA= -40° to PA= -65° 
with increasing distance from the nucleus. 



sorbers in QSOs peaks around the UV emission-line velocity. 
Nevertheless, one might expect to find gas with random motions, 
both blue- and red-shifted, relative to the host galaxy intercept- 



ing our line of sight either within the host galaxy itself or asso- 
ciated with its neighbors. Indeed, one component of the absorp- 
tion in the Seyf ert 1 galaxy Mrk 279 has a positive velocity of 
+90 km s -1 , and lScott et al.l (12004) have suggested that it is inter- 
stellar material in the host galaxy, or gas from a recent encounter 
with a close companion, perhaps similar to a high-velocity cloud 
(HVC) in our own Milky Way. 

Since our absorption-line observations only probe gas that 
is along the line of sight, taking a broader view of the en- 
vironment of the active nucleus in Mrk 509 should be help- 
ful. The optical imaging and two-dimensional optical spec- 
troscopy of Phillips et al.l d!983l) reveal a compact host, perhaps 
of type SO. The spectroscopy shows kinematically distinct low- 
ionization and high-ionization emission-line components. The 
low-ionization gas covers the face of the galaxy and rotates as 
a disk about an axis at a position angle (PA, measured East of 
North) of -45°. The high-ionization gas also extends to more 
than 10" from the nucleus, and it is mostly blueshifted, with 
velocities as high as -200 kms~', although there are also re- 
gions SE of the nucleus with positive velocities. Radio images 
of Mrk 509 show elongation and ext ensions to both the NE and 
the SW dSingh & Westergaardl 19921) . 
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To illustrate these features on the spatial scale of our HST 
observations, we have retrieved the WFPC2 F547M image 
of Mrk 509 obtained by iBentz et all (120091) from the MAST 
archive. We used MultiDrizzle dFruchter et al.ll2009l) to combine 
the longest PC exposures (with exposure times of 60 s, 160 s, 
and 300 s) and eliminate cosmic rays. Figure [9] shows the re- 
sultant image (covering the inner 4 kpc of Mrk 509). Note that 
this image does not capture the full extent of Mrk 509. The 
2.5"-diameter COS aperture covers only the central 870 pc of 
Mrk 509, while the most prominent features in the HST im- 
age, possibly an irregular starburst ring, have an outer radius of 
< 2 kpc. If this ring-like feature is due to recent star formation, 
this may be the source of ionizati on for the low-ioni zation ro- 
tating disk in the spectroscopy of Phill ips et al.l dl983l) . There is 
also an unusual, linear feature to the NE, which is not aligned 
with the nucleus, although it lies in the general direction of the 
elongations to the NW in the radio maps. Deep ground-based 
CCD images of Mrk 5 09 show an amorphous optical morphol- 
ogy dMacKentvl 1 1 990h with a broad extension to the SE in the 
same direction as the elongation in the 6-cm radio image of 
Singh & Westergaard ( 1992). Although the light is highly con- 
centrated, MacKentyJ (119901) traces an exponential disk out to a 
radius of more than 30 kpc. 

In addition to this macroscopic view of the nuclear region, 
optical spectropolarimetric observations have probed structures 
that are unresolved in the optical or radio images. Optical spec- 
tropolarimetry of Mrk 509 suggests that there are multiple scat- 
tering regions, each consistent with various aspects of the unified 
model of AGN (Young 2000). The optical continuum is polar- 
ized at PA=139-154°, along the same direct ion as the inner por- 
tions of the 6-cm and 20-cm radio images (Singh & Westergaard 



1992). The po larization is also variable (Young et al. 1999), and 
Young (2000) identified this as radiation scattered from the ac- 
cretion disk. The Ha emission line exhibits more complexity in 
its polarization. The core is polarized at PA=60°, nearly perpen- 
dicular to the radio jet, as observed in Type II AGN, where the 
scattering surface is the mirror of hot electrons in the cone along 
the radio axis, perpendicular to the plane of the obscuring torus. 
The wings of Ha are polarized at PA=168°, similar to, but not 
quite aligned with the inner radio jet. Figure [9] shows these ori- 
entations relative to other structures in Mrk 509. 

As we noted in §3.4, the COS and FUSE spectra show 
significantly more residual light at the bottoms of the intrin- 
sic absorption-line troughs than does the STIS spectrum. Since 
COS and FUSE used larger apertures, we argued that much 
of this light was scattered light, primarily from a region less 
than 0.44445" in diameter surrounding the active nucleus. (This 
scattered light is from a much smaller region than the broad 
Balmer emission that Mediavilla et al.l (Il998l) detected using 
two-dimensional spectroscopy over a region of ~ 8" surround- 
ing the nucleus, which they interpret as scattered light from the 
AGN.) A distance of 0.225" corresponds to a radial distance of 
< 150 pc, which would place it in the typical region associated 
with reflective hot-electron mirrors in Type II AGN, as we noted 
in §3.4. This is consistent with the polarization component seen 
in the core of the Ha emission line that is polarized perpendicu- 
larly to the radio axis. The amount of scattered light we see, ap- 
proximately 5% of the total flux, is also typical of the amount of 
scattered light seen in T ype II AGN like NGC 1068 dMiller et all 
ll99lHAntonucrill993l) . 

Thus our line of sight to the nucleus in Mrk 509 may pass 
through a complicated mix of components. Based on the kine- 
matics and physical characteristics of the absorption components 



revealed by our spectrum, we will try to link the absorption com- 
ponents to the various features in Mrk 509. 

4.2. Characterizing the Absorbers 

In this paper we take a kinematic approach to identifying 
and characterizing the absorbers in Mrk 509 (as oppose d to 
the ionization-stat e appr oach used by IDetmers et al.ll201ll and 
Ebre ro et al.ll2.01 lb . As in lEbrero et al.l ( 1201 lb . our discussion is 
a preliminary assessment of the properties of the absorbers. A 
more comprehensive analysis of the X-ray and UV properties of 
the absorbing gas is deferred to a future paper (Ebrero, in preps.). 

On the largest velocity scales, the intrinsic absorption in 
Mrk 509 can be characterized by two main troughs. The 
most blue-shifted trough, extending from -200 kms~' to 
-450 kms~', corresponds to the velocities of the most promi- 
nent blue-shift ed absorption comp onent in the RGS spectrum 
(-319 kms~'; Detmers et al. 2011) and the Chandra spectrum 
(-196 km s _1 ; Ebrero et al. 201 1). (This velocity offset between 
the RGS and LETGS results is not significant, and it lies within 
the range of systematic errors in the relati ve RGS and LE TGS 
wavelength scales as discussed by Kaas tra et al.l 1201 lal and 
Ebre ro et al.ll201 ll ) Given the variability we have seen in por- 
tions of this trough and its blueshift, we identify this trough as 
an outflow from the AGN. 

The other main absorption trough in Mrk 509 extends from 
-100 kms -1 to +250 kms -1 . This also matches the kinemat- 
ics of the other main absorption component detected in the 
X-ray spectra, with velocities meas ured in the RGS spectrum 
of -13 kms 1 (IDetmers et al.ll201 ll) a nd of +73 kms~' in the 
Chandra spectrum (Ebr ero et al.ll20 1 1). The ionization level in 
the X-ray for this component is lower, and this gas is not in 
pressur e equilibrium with t he gas in the blue-shifted absorption 
trough dEbrero et al.ll201 11) . As we discuss below, the compo- 
nents in this trough may arise in several distinctly different re- 
gions of Mrk 509. The higher spectral resolution of our UV spec- 
trum permits us to dissect these absorption features in greater de- 
tail, so we now discuss the properties of individual components 
that we have identified in the UV absorption spectrum. 



4.2.1. Components 1a, 1, and 1b 

Component la is detected only in the far wing of the Lya absorp- 
tion trough. There is a distinct depression here that is not seen in 
the profiles of any of the other absorption lines, including O vi. 
This suggests that the gas here may be very highly ionized, per- 
haps representing the only indication in the UV of a more exten- 
sive outflow at higher velocities that is see n in some of the most 
highly ionized X- ray absorption features (IDetmers et al.ll201lt 
lEbrero et al.ll20H . 

Components 1 and lb show definite evidence of variability 
from comparing the STIS and COS observations. Both vary in 
Lya, and component 1 varies in N v as we discussed in §3.4. 
Although more closely spaced in time, the FUSE observations 
do not permit us to set any interesting limits on the density or 
location of the gas in these components. These components may 
have varied between the two FUSE observations, but not at a 
level of significance that permits us to infer densities or dis- 
tances. The variations in N v and Lya are consistent with ion- 
ization changes in response to a decrease in the continuum flux. 
Recombination times for hydrogen are too slow to permit us 
to set any interesting limits on the density. However, N v re- 
combines nearly two orders of magnitude faster, and it provides 
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greater lev erage in sensing the de nsity of the absorbing me dium. 
Following iKrolik & KrissT (1 19951) and iNicastro et all i 19991) . re- 
combination and ionization time scales depend not only on the 
density, but also on the relative populations of the ionization 
states involved: t rec = (ni/n; + i)/(n e a rec ). Using the photoion- 
ization model for the Mrk 509 absorbers from iKraemer et al.l 
d2003l) . component 1 has an ionization parameter log ^ = 0.67 
(where £ = L; on /(nr 2 )). At this level, the ionization fraction for 
N v is ~ 0.4, and the relative equality of populations among 
the neighboring ionization states of nitrog en makes ioniz ation 
and recombination timescales similar (Nicastro et al. 1999). For 



4.2.4. Component 4a 



a rec - 8.96 x 10 cm i s~ (Nahar 2006) at a temperature 
of 20,000 K, and a time between the STIS and COS observa- 
tions of 2.733 x 10 8 s, we get a lower limit on the density of 
n e > 160 cm' 3 . Sin ce we are using the photoionization mod- 
els of Kraemer et al. (2003), we use an ionizing luminosity from 
their SEP in lYaqoob et all d2003l) . L ion = 7.0 x 10 44 ergs s _1 . 
Together with the ionization parameter log£ = 0.67, this gives 
an upper limit on the distance of r < 250 pc. 

Component 1 is well defined in both N v and C iv, but in 
Lya and O vi it is merely part of a smoothly varying trough that 
runs from the deepest point i n components 2 and 2b. As mod- 
eled by IKraemer et all (120031) (who only used the STIS data), 
this trough has only moderate ionization, probably even lower 
than what they model since they grossly overpredict the amount 
of O vi that should be present at this velocity. (Their trough at 
-420 km s is three times deeper than the observed O vi /U032 
absorption in the FUSE spectrum at this velocity.) In fact, as we 
discussed in §3.4, given the relative character of the variabil- 
ity in Lya and N v, the ionization state of this feature is below 
the peak in ionization for N v. Component lb appears to have 
higher ionization since it is strong in O vi but very indistinct in 
N v. Although moderately strong in Lya, it is also weaker than 
component 1 . If either of these components are associated with 
the X-ray portion of the outflow (which they do not match in 
detail kinematically), they are likely to be higher-density, lower- 
ionization clumps embedded in a more tenuous, highly ionized 
outflow. 



4.2.2. Components 2 and 2b 

Components 2 and 2b represent the deepest portions of the 
blue half of the absorption trough in the Mrk 509 spectrum. 
Kinematically, these components are the closest match to one 
of the major components detected in the X-ray spectrum of 
Mrk 509. In the XMM -Newton/RGS spectrum, component 2 of 
Det mers et alJ (1201 lh is almost an exact match in velocity to UV 
component 2 (-319 kms' 1 vs. -321 kms' 1 , respectively). The 
ionization state of the gas detected in the UV absorption lines, 
however, is not particularly high, and its total column density is 
less than that seen in the X-ray. Again, this may be an example of 
a higher-density, lower-ionization clump embedded in the more 
highly ionized outflow seen in the X-ray. 



4.2.3. Component 3 

Component 3 is the closest match i n velocity to Comp onent 2 
of the Chandra LETGS spectrum dEbrero et alJ 1201 II) . but in 
the UV it too has rather low ionization, similar in character to 
UV component 1 . As with the previous components, this gas de- 
tected in the UV is likely to be a higher-density clump embedded 
in more highly ionized outflowing gas. 



Component 4a, first identified by Kra emer et al.l (120031) . shows 
up most clearly in N v and C iv. In velocity it is close to compo- 
nent 1 identified in the RGS spectrum by Detmers et al. (201 1). 
Its modest outflow velocity could suggest an association with 
the base of the outflowing wind in Mrk 509. However, like all 
the components above, it is rather low in ionization, and its col- 
umn density is lower than that seen in the X-ray absorber. It is 
probably not directly associated with the X-ray absorbing gas. 

4.2.5. Component 4 

Perhaps the easiest component to characterize is #4. This ab- 
sorber has the lowest ionization state of any of the components. It 
is the only one in which Si iv is detected, and it has th e strongest 
Cm absorption in the original FUSE spectrum of iKriss et al.l 

(2000). Its velocity of 22 km s~ is the closest to our adopted 

systemic redshift for Mrk 509. Given the high column densi- 
ties, relatively low ionization state and velocity coincidence, it 
is likely that this co mponent is the interst ellar medium of the 
Mrk 509 host galaxy. Philli ps et alJ (1 19831) find a systemic ve- 
locity for the rotating low-ionization gas disk in their observa- 
tions that is in good agreement with the velocity of the [O in] 
line-center velocity which they adopted as the systemic velocity 
of the host galaxy. To within their errors of ±30 kms' 1 , all of 
these velociti es are consistent with the velocity we measure for 
component 4. IKraemer et alJ d2003l) find that the broad range of 
UV ionization states present at this velocity requires more than 
one component in a photoionization model solution. Component 
4 has the closest match in velocity to compon ent 1 identified 
in the RGS spectrum by iDetmers etaD (1201 ll) . Within the er- 
rors of the Chandra spectrum, Component 4 may also be associ- 
ated with component 1 from the LETGS spectrum ( Ebrero et al. 
1201 ll) . Given its low ionization parameter, Ebr ero et al.l (1201 lh 
suggest that this X-ray absorption may be associated with the 
ISM or halo of the host galaxy of Mrk 509. 



4.2.6. Component 5 

Component 5 is similar in character to Component 4, but it has 
lower overall column density and slightly higher ionization. Its 
slightly positive velocity is a good match to component 1 in the 
LETGS spectrum. Highly ionized gas at this velocity is also seen 
in the three-dimensional spectroscopy of Phillips et al. (1983), 
and it is within the velocity profile of the echelle spectrum of 
[O m] in the nuclear region of Mrk 509. It appears to have varied 
in strength in H i between the two FUSE observations. As it is 
difficult to disentangle this component's variations from those of 
components 6a, 6, and 6b, we discuss the implications of this 
variation in the next subsection. 



4.2.7. Components 6a, 6, and 6b 

Although component 6 has a positive velocity, this does not pre- 
clude it from an association with the active nucleus in Mrk 509. 
In fact, if one examines the [O in] e chelle spectrum of th e nu- 
clear region in Mrk 509 in Fig. 3 of lPhillips et alJ d 19831) . one 
can see that the actual peak of the line is at a positive velocity 
of +120 kms' 1 , matching our velocity for absorption in com- 
pone nt 6. The kinematics of X-ray-heated winds modeled by 
Bals ara & Krol ik ( 1993) shows that portions of the flow that are 
evaporated off the inner edge of the obscuring torus can be cap- 
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tured by the central black hole. This results in streamlines with 
positive, inflow velocities. 

The variability in H i absorption we see over the 10 months 
between the two FUSE observations suggests a close relation- 
ship to the active nucleus for components 5 and 6, and, as for 
Component 1 , it allows us to set an upper limit on the distance 
of this gas from the central engine. From 1999 November to 
2000 September, over an interval of 303 days, the continuum 
flux from Mrk 509 dropped by 30%, and the hydrogen column 
density in components 5 and 6 increased by 14%, as we showed 
in §3.4. This gives us an upper limit on the recombination time 
for the gas, and hence a lower li mit on the dens i ty. A gain using 
the photoionization model from iKraemer et all (2003), compo- 
nent 6 has an ionization parameter log £, - 0.71, and an ion- 
ization fraction for hydrogen of 2.3 x 10~ 5 . For a time interval 
between the two FUS E observations of 2.62 X 10 7 s, and a rec = 
1.43 x 10~ 13 cm 3 s _1 dOsterbrock & Ferland 2006) at a tempera- 
ture of 20,000 K, we get a limit on the density of n e > 6.0 cm -3 . 
The ionizing luminosity of Li on = 7.0 x 10 44 ergs s -1 together 
with the ionization parameter log£ = 0.71, gives an upper limit 
on the distance of r < 1.5 kpc. So, component 6 is definitely 
within the confines of the host galaxy, and likely in the near 
vicinity of the nucleus. 

4.2.8. Components 7a and 7 

Component 7 at +219 kms~', from its velocity, is not outflow- 
ing, but infalling to Mrk 509. Component 7a is strongly blended 
with #7 in O vi and Lya; it only shows up as a kinematically 
independent feature in N v. Neither component shows any sig- 
nificant variability either between the two FUSE observations 
in 1999 and 2000, or between the longer interval between STIS 
and COS from 2001 to 2009. Given that it is optically thin in 
N v, C iv, and Lya, this lack of variability is not being masked 
by saturation effects that limit our ability to see variations in 
other components. Based on just the qualitative appearance of 
the UV spectra, this is the highest ionization parameter trough in 
the whole spectrum — O vi is very strong, N v is present, and C iv 
is not detectable. The stacked RGS spectrum at these velocities 
(Detmers et al. 201 1) does not reveal absorption in O vi Ka, al- 
though the upper limits are consistent with the FUSE spectrum. 
So, even though the ionization level exhibited in the UV is high, 
it is st ill not as high a s the blue-sh ifted trough in the X-ray spec- 
trum. IKraemer et all d2003l) and iKriss et all d2OOO0 found their 
highest ionization parameters for this component as well. 

We suggest that components 7 and 7a might be similar to 
HVC complexes seen in our own galaxy, except that the ion- 
ization of these clouds in Mrk 509 would be completely domi- 
nated by the radiation from the active nucleus rather than the ex- 
tragalactic background radiation and starlight from the galaxy. 
These clouds may be infalling extragalactic material, or they 
may be tidal material stripped from an encounter with a com- 
panion in the past. Although Mrk 509 does not have an obvious 
intera cting companion, its disturbed outer isophotes (MacKentv 



1990) suggest there may have been an encounter in the more 
distant past. iThom et al] J2008I) find that HVC Complex C in 

the Milky Way has an average density log n 2.5, dimensions 

3x15 kpc, a distance of 10 + 2.5 kpc, and a total mass of 8.2+ 4 gX 
10 6 M Q , If components 7 and 7a in Mrk 509 have roughly similar 
densities of n - 0.01 cm" 3 , the n we can use the photoionization 
results of IKraemer et al J (2003) to esti mate their s i ze and dis- 
tance from the nucleus of Mrk 509. IKraemer et al. (2003) give 
a total column density of Nh — 4 x 10 19 cirT 2 for component 



7, which implies a size (as a cube) of 1.3(0.01 cirT 3 /n) kpc. For 
their quoted ionization parameter of log £, - 1.33, and an ioniz- 
ing luminosity of L; on = 7.0 x 10 44 ergs s , the distance from the 
center of Mrk 509 is 19 (0.01 cirT 3 /n)° 5 kpc. This is just about 
right for a structure comparable to a Milky Way HVC. For a den- 
sity this low and ionization this high, the recombination time for 
H i in components 7a and 7 would be 80(0.01 cirT 3 /n) years, so 
the lack of variability in our observations is consistent with this 
hypothesis. 

4.2.9. Smoothly distributed gas in the outflow 

While we have characterized the absorption features in Mrk 509 
in terms of discrete kinematic components, the smoothly vary- 
ing profiles of the Lya and O vi absorption in the blue-shifted 
outflow trough (see Fig. [5]) suggest that additional gas may be 
present that is smoothly distributed over a range in velocities. 
The tail to higher velocities observed as component la in Lya 
may be a part of this distribution. If this smoothly distributed 
gas is very highly ionized, one might expect that the only traces 
remaining in the UV absorption spectrum would be the most 
highly ionized species, like O vi, or the most abundant ones, like 
H i. This more smoothly distributed, highly ionized gas could be 
the counterpart of the features detected in the RGS and LETGS 
X-ray spectra. 

The discrete UV absorption components we have identified 
show up most clearly in N v and C iv. As noted in §3.3, these of- 
ten do not precisely align with similar features in O vi and Lya. 
For example, in the region surrounding component lb, there is 
little C iv and N v absorption, but moderate depth in O vi and 
the hydrogen lines, although it is a local minimum in H i. The 
lower ionization, more discrete components we observe in N v 
and C iv could be dense, lower-ionization clumps embedded in 
this flow, or clouds in the interstellar medium of Mrk 509 that ar e 
being entrained in the flow. As shown by Ebre ro et al.l d2011l) . 
these lower-ionization components are not in pressure equilib- 
rium with the other X-ray absorbing components. Such discrete 
clouds with similar kinematics are spatial ly resolved and ob- 
served in the nearest A GN, e . g., NG C 4151 dDas et alj2005l) and 
NGC 1068 dDas et alJ l2006. 2007). This outflowi ng gas is also 
visible in the three dimensional spectroscopy of Phillips et al. 
(1983). 



4.3. Thermal Winds vs. Disk Winds 

Disk winds are a popular explanation for outflows from AGN 
(Konigl&Kartie 1994; Murrav & Chiang 1995; Elvis 200'~ 
Progaetal. 2000; Richards et al. 2010; Fukumura et al. 201 
For the high velocities observed in broad-absorption-line (BAL) 
QSOs, an origin for the wind from the accretion disk deep in 
the potential well of the central black hole is a natural infer- 
ence. Fo r the low velocities w e observe in Mrk 509, it is not 
so clear. Krae mer et al.l (120031) summarizes the problems faced 
in trying to explain the outflow in Mrk 509 with a disk wind. 
Briefly, Mrk 509 is likely to be observed nearly face on to the 
accretion disk based on the polarized flux in the broad wings of 
Ha that is polarized parallel to the radio jet, suggest ing an ori- 
gin in scattering from the face of the accretion disk (1 Young et al.l 
1999; Young 2000). Since disk winds roughly lie in the plane of 
the disk, the flow would be predominantly perpendicular to our 
line of sight. This can explain the low velocities since most of 
the motion would be transverse to our line of sight, but it makes 
it difficult to explain the near-unity covering factors of the ab- 
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sorbers if the UV continuum is produced interior to the launch 
point of the wind. 

The modest velocities of hundreds of kms -1 we observe 
here are more characteristic of thermal winds produced in the 
AGN environment. Irradiation of the obscuring torus in an 
AGN by the nuclear flux produces material that fills the cone 
above the torus and flows outward at velocities of hundreds of 
kms , typical of the sound speed in the gas (Balsara & Krolik 
119931: lKrohk&Krisslll995L 1200 lb . The presence of some red- 
shifted components is also consistent with portions of an X- 
ray-heated wind that a re captured by the central black hole 
(Balsara & Kjolikf l993). These winds are also expected to have 
multiple phases of gas as the surface o f the torus is ablated 
by th e radiation from the central engine dKrolik & Krissl [19951. 
2001). The most highly ionized portions of this flow can pro- 
duce the X-ray absorption. The dense part of the wind near the 
base of the ionization cone may be the slightly extended scat- 
tering region producing the extra flux we see in the bottoms of 
the absorption troughs and also produce the polarization seen 
in the core of the Ha emission line. As we noted in the prior 
section, the lower-ionization UV absorption components could 
be denser clumps embedded in this flow. The extended high- 
ionization [Oii i] emission outflow ing from the nuclear region 
and studied by Phillips et al. (1983) could be analogous to the 
discrete clouds imaged with HST in nearby AGN that show sim- 
ilar kinematics. 



5. Conclusions 

We have presented HST/COS observations of the Seyfert 
1 galaxy Mrk 509 obt ained as part of an extensive mul- 
tiwavelength campaign dKaastra et al.l 1201 lbl) on this galaxy 
in late 2009. Our UV spectra were obtained simultaneously 
with Chandra LETGS spectra (Ebreroetal. 2011). The cam- 
paign also included immediately prior observations with XMM- 
Newton, Swift, and INTEGRAL. Our spectra cover the 1155- 
1760 A wavelength range at a resolution of ~ 15 kms 1 , and 
they are the highest signal-to-noise observations to date of the 
intrinsic absorption components of Mrk 509. This enables us to 
trace additional complexity in the absorption tro ughs compared 
to prior STIS observations dKraemer et al.l l2003). As part of our 
analysis, we have also examined archival FUSE observations, 
which provide additional information on variability characteris- 
tics of the absorbers. 

In our COS spectra we identify a total of 14 kinematic com- 
ponents. Six of these are blends that represent subcomponents 
of previously known features. The most blueshifted portions of 
the Lyo- absorption trough (components 1 and lb) show a signif- 
icant decrease in column density compared to prior STIS ob- 
servations. At the same time, component 1 in N v increased 
in strength compared to the STIS observations. Using recom- 
bination timescale arguments, the variation in the N v column 
density allows us to set a lower limit on the gas density of 
> 160 cm 4 , and an upper limit on the distance of < 250 pc 
for this most blue-shifted component. We also detect variabil- 
ity between the two FUSE observations from 1999 November 
(Kriss et al. 2000) to 2000 September in component 6 at a pos- 
itive velocity of +120 kms 1 , which corresponds to the peak 
of the [O ml emissio n-line profile from the echelle spectrum of 
Phillips et al. ( 1983). Based on the recombination time, we set a 
limit on the density of this component of > 6 cm 3 , and set an 
upper limit on its distance from the nucleus at < 1 .5 kpc. 

The COS spectra also show greater residual light at the bot- 
toms of the absorption troughs relative to the STIS spectrum. 



The fraction of residual light seen in COS is comparable to 
that seen in the FUSE spectra. We attribute this residual light 
to scattering from an extended region near the nucleus, which 
is consistent with the larger apertures used for the COS and 
FUSE observations relative to the STIS spectrum. Analysis of 
the cross-dispersion profile in the COS spectrum limits the spa- 
tial extent of this scattering region to < 0.5", or < 170 pc. 
Scattering from such a region would be consistent with electron 
scattering from the base of an ionization cone in Mr k 509, and it 
would explain the polarization in the core of Ho- ( Youn get al.l 
1999) which is observed to be perpendicular to the radio jet 
(Sing h & Westergaardlll992b . 

The velocities of the two main absorption troughs roughly 
correspond to the kinematics of the detected X-ray absorption. 
The UV absorption in Mrk 509 arises from a variety of sources. 
The lowest-velocity absorption trough covers a velocity range of 
- 1 00 km s~ ' to +250 km s~ ' . The deepest portion of this trough 
in the UV is at the systemic velocity of Mrk 509, and it has the 
lowest ionization. This is also true for the X-ray absorption. We 
attribute this portion of the absorption to the interstellar medium 
or galactic halo of Mrk 509. The most redshifted portion of the 
trough has characteristics comparable t o the high- velocity cloud 
Complex C in our own Milky Way dThom et al.l2 008). At a den- 
sity of 0.01 cm 4 (which we can only assume, not measure), it 
would have a size of roughly 1.3 kpc at a distance of 19 kpc from 
the center of Mrk 509. 

The most blue-shifted absorption trough in Mrk 509 is an 
outflow from the active nucleus. The velocities in this trough 
of -200 kms ' to -450 kms 1 correspon d to the most highly 
ionized portion of t he X-ray absorbing gas (Det mers et alj2011t 
lEbrero et al.l 1201 lb . and they overlap with the extended emis- 
sion from high-ionization gas that cove rs the face of Mrk 5 09 in 
the three-dimensional spectroscopy of Philli ps et al.l (11983). The 
outflow velocities, embedded clumps of lower ionization gas, the 
presence of an extended scattering region, an d the extended out- 
flowing emission seen in [Oiii] (Phillip s et all 19831) are all con- 
sistent with an origin in a multiphase thermal wind produced by 
the irradiation of an obscuring torus by the active nucleus. 
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Feature Aq Flux 
(A) (Kr' 4 ergcnr 2 s-' A' 1 ) 



v sys 
(kms- 1 ) 



FWHM 

(kms" 1 ) 



Lya 


1215.670 


Lya 


1215.670 


Lya 


1215.670 


Lyce 


1215.670 


Nv 


1238.821 


Nv 


1242.804 


Sin 


1260.420 


O i+Si ii 


1303.57 


Cn 


1335.205 


Si iv 


1393.755 


Si iv 


1402.770 


Civ 


1548.195 


Civ 


1550.770 


Civ 


1548.195 


Civ 


1550.770 


Civ 


1548.195 


Civ 


1550.770 


Civ 


1549.050 


Niv] 


1486.496 


Sin 


1526.710 


Hen 


1640.480 


Hen 


1640.480 


Om] 


1663.445 



9.1 ±2.5 
73.9 ± 2.5 

410.0 ±1.8 
838.0 ±1.8 
22.1 ±0.6 
11.1 ±0.2 
65.3 ±1.0 
31.5 ±0.6 

8.0 ±1.5 
85.0 ±0.2 
42.5 ±0.1 

4.2 ±1.4 

2.1 ±0.3 
78.5 ±0.5 
39.3 ±0.5 
199.0 ±1.9 

99.3 ±1.9 
461.0 ±1.1 

7.0 ±0.1 
15.7 ±1.3 
8.5 ±0.3 

55.4 ±0.5 
66.0 ±1.6 



166 ±31 


300 ± 23 


-151 ± 15 


1339 ± 33 


-235 ± 8 


3039 ± 14 


-156 ±7 


9677 ± 26 


205 ± 22 


2164 ± 45 


205 ± 22 


2164 ±45 


-707 ± 36 


5652 ± 66 


36 ±33 


3307 ± 90 


8 ±67 


2176 ± 424 


492 ±9 


5084 ± 36 


491 ±9 


5084 ± 36 


-204 ± 39 


300 ± 39 


-204 ± 39 


300 ± 39 


-206 ± 7 


1974 ± 15 


-206 ± 7 


1974 ± 15 


-206 ± 7 


4233 ± 15 


-206 ± 7 


4233 ± 15 


158 ±10 


10299 ± 56 


-33 ±61 


2181 ±108 


2 ±27 


2181 ±108 


-113 ±39 


1153 ±83 


-199 ±25 


4065 ± 61 


83 ±19 


3900 ± 69 



Table 4. Emission Features in the STIS Spectrum of Mrk 509 



Feature 



(A) 



Flux 

(lO^ergcm^s- 1 A" 1 ) 



v sys 
(kms- 1 ) 



FWHM 

(kms- 1 ) 



Lya 


1215.6700 


Lya 


1215.6700 


Lya 


1215.6700 


Lyce 


1215.6700 


Nv 


1238.8210 


Nv 


1242.8040 


Sin 


1260.4200 


O i+Si ii 


1303.5700 


Cn 


1335.2050 


Si iv 


1393.7550 


Si iv 


1402.7700 


Civ 


1548.1950 


Civ 


1550.7700 


Civ 


1548.1950 


Civ 


1550.7700 


Civ 


1548.1950 


Civ 


1550.7700 


Civ 


1549.0530 


Niv] 


1486.4960 


Sin 


1526.7070 


Hen 


1640.4800 


Hen 


1640.4800 



6.8 ±5.5 
102.0 ± 0.9 
350.0 ± 2.9 
531.0 ±4.4 
17.4 ±0.3 
8.7 ± 0.2 
36.1 ±1.0 

21.8 ±2.5 
4.7 ± 2.4 

54.9 ± 3.0 

27.4 ±1.5 
6.7 ±1.8 
3.4 ±0.9 

65.9 ± 3.6 

33.0 ±1.8 
148.0 ±2.0 

74.1 ± 1.0 
389.0 ±1.5 

12.3 ± 0.6 
15.8 ± 1.4 
5.4 ±0.8 

39.5 ±3.1 



166 ±575 
85 ± 15 
-290 ±4 
-11± 19 
371 ±44 
371 ±44 
-921 ± 130 
189 ±62 
106 ± 284 

381 ± 194 

382 ±194 
-204 ±151 
-204 ±151 

-103 ±7 
-103 ±7 
-103 ±7 
-103 ±7 

-605 ± 10 
262 ±71 

-175 ±86 
38 ±48 

-636 ± 67 



300 ±3 
1712 ±66 
3483 ± 25 
10752 ± 56 
2089 ± 39 
2089 ± 39 
5928 ±151 
3115 ±94 
1555 ± 459 
4374 ± 97 
4374 ± 97 
300 ± 95 
300 ± 95 
2108 ± 30 
2108 ±30 
4102 ±93 
4102 ±93 
9392 ± 156 
3795 ± 20 
3795 ± 20 
888 ± 135 
4824 ± 104 
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Table 5. Emission Features in the 1999 FUSE Spectrum of Mrk 509 



Feature 


^0 


Flux 


v sys 


FWHM 




(A) 


(lO^ergcm^s- 1 A" 1 ) 


(kms- 1 ) 


(kms- 1 ) 


Cm 


977.0200 


60.0 ±4.7 


-8 ±25 


4000 ± 172 


Nm 


989.7990 


56.1 ±1.8 


-546 ±51 


8917 ±217 


Ly/J 


1025.7220 


0.6 ±5.1 


-812 ± 197 


8917 ±217 


Ovi 


1031.9260 


100.0 ±3.2 


-812 ± 197 


8917 ±217 


Ovi 


1037.6170 


50.1 ±1.6 


-812 ± 197 


8917 ±217 


Ovi 


1031.9260 


69.0 ±3.1 


193 ± 48 


4182 ± 99 


Ovi 


1037.6170 


34.5 ± 1.6 


193 ± 48 


4182 ± 99 


Ovi 


1031.9260 


3.4 ±0.6 


-203 ± 33 


300 ± 47 


Ovi 


1037.6170 


1.7 ±0.3 


-203 ± 33 


300 ± 47 



Table 6. Emission Features in the 2000 FUSE Spectrum of Mrk 509 

Feature T Flux v^ FWHM 

(A) (IP' 14 erg cm' 2 s" 1 A' 1 ) (km s"') (km s"') 

858 ± 665 9157 ± 860 

137 ±105 729 ±277 

-297 ±245 7844 ±187 

-566 ±245 7844 ±187 

-565 ±245 7844 ±187 

-564 ±245 7844 ±187 

199 ± 855 2585 ± 272 

199 ± 855 2585 ± 272 

-203 ± 129 300 ± 250 

-203 ± 129 300 ± 250 



Cm 


977.0200 


Cm 


977.0200 


Nm 


989.7990 


Ly/J 


1025.7220 


Ovi 


1031.9260 


Ovi 


1037.6170 


Ovi 


1031.9260 


Ovi 


1037.6170 


Ovi 


1031.9260 


Ovi 


1037.6170 



26.8 ± 


1.1 


2.5 ± 


0.1 


23.0 ± 


1.6 


0.1 ± 


1.2 


109.0 ±5.0 


54.3 ± 


:2.5 


43.2 ± 


4.0 


21.6 + 


2.0 


4.2 ± 


1.0 


2.1 ± 


0.5 
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Table 9. Intrinsic Absorption Features in the COS Spectrum of Mrk 509 



Ion 


Component 


N 


v sys 


b 


fc 




# 


(10 12 cm" 2 ) 


(km s" 1 ) 


(km s" 1 ) 




Hi 


la 


6.3 ± 0.7 


-426 ± 5 


13 ±1 


0.802 ± 0.087 


Hi 


1 


49.8 ±4.9 


-385 ± 7 


30 ± 1 


0.802 ± 0.087 


Hi 


lb 


2.2 ±2.1 


-361 ± 8 


12 ±2 


0.802 ± 0.087 


Hi 


2 


349.7 ± 87.5 


-316 ±6 


25 ± 1 


0.965 ± 0.084 


Hi 


2b 


133.2 ±66.1 


-283 ± 5 


26 ±3 


0.965 ± 0.084 


Hi 


3 


75.6 ± 9.2 


-240 ± 7 


32 ± 1 


0.843 ±0.126 


Hi 


4a 


262.0 ±9.1 


-59 ±5 


29 ±1 


0.953 ± 0.005 


Hi 


4 


134.9 ± 27.4 


-8 ±6 


31 ±1 


0.984 ± 0.028 


Hi 


5 


217.3 ± 30.0 


44 ±8 


42 ±5 


0.970 ± 0.024 


Hi 


6a 


34.6 ± 10.0 


132 ±7 


33 ±2 


0.993 ± 0.024 


Hi 


6 


312.5 ±35.6 


125 ±6 


33 ±2 


0.993 ± 0.024 


Hi 


6b 


0.1 ± 1.0 


166 ±8 


33 ±2 


0.993 ± 0.024 


Hi 


7a 


53.8 ± 12.8 


185 ±5 


16 ±3 


0.861 ±0.051 


Hi 


7 


123.9 ±2.5 


220 ±5 


32 ±2 


0.861 ±0.051 


Nv 


1 


107.0 ±9.5 


-408 ± 5 


10 ±1 


0.712 ±0.026 


Nv 


lb 


17.5 ± 1.6 


-375 ± 7 


21 ±3 


0.678 ± 0.010 


Nv 


2 


149.0 ± 10.7 


-321 ±5 


13 ±1 


0.950 ± 0.017 


Nv 


2b 


130.6 ±11.3 


-291 ± 6 


19 ±2 


0.729 ± 0.025 


Nv 


3 


88.6 ±8.1 


-244 ± 5 


16 ±1 


0.698 ± 0.045 


Nv 


4a 


93.6 ± 6.8 


-59 ±5 


13 ±1 


0.945 ±0.031 


Nv 


4 


323.6 ± 12.2 


-19 + 5 


14 ±1 


1.000 ±0.012 


Nv 


5 


38.3 ±11.1 


37 ±5 


13 ±1 


1.000 ±0.250 


Nv 


6a 


20.6 ± 6.4 


79 ±7 


41 ±12 


1.000 ±0.023 


Nv 


6 


56.1 ±4.0 


121 ±5 


13 ±1 


1.000 ±0.023 


Nv 


6b 


5.5 ± 1.8 


147 ±8 


11 ± 3 


1.000 ±0.023 


Nv 


7 


21.0 ±2.9 


184 ±6 


16 ±1 


0.836 ± 0.094 


Nv 


7a 


23.5 ± 2.7 


222 ±6 


16 ±1 


0.836 ± 0.094 


Si iv 


4 


2.5 ± 0.2 


-22 ±6 


10 ±1 


1.000 ±0.146 


Civ 


1 


31.2 ±1.5 


-406 ± 6 


14 ±1 


0.892 ± 0.026 


Civ 


lb 


16.3 ±6.1 


-365 ± 5 


14 ±1 


0.935 ±0.261 


Civ 


2 


136.3 ±41.7 


-320 ± 6 


13 ±1 


0.979 ± 0.035 


Civ 


2b 


128.9 ±41.4 


-292 ± 7 


22 ±3 


0.896 ± 0.030 


Civ 


3 


44.7 ± 7.2 


-245 ± 5 


22 ±2 


0.987 ± 0.023 


Civ 


4a 


66.3 ±2.1 


-60 ±5 


12 ± 1 


0.961 ±0.031 


Civ 


4 


250.0 ±11.0 


-18 ±5 


16 ±1 


0.967 ± 0.007 


Civ 


5 


36.2 ± 16.9 


34 ±5 


19 ± 1 


0.939 ±0.109 


Civ 


6a 


5.6 ±2.3 


88 ±7 


12 ±2 


0.986 ± 0.049 


Civ 


6 


12.3 ± 6.4 


122 ±5 


10 ±1 


0.986 ± 0.049 


Civ 


6b 


17.6 ±5.3 


120 ±8 


24 ±4 


0.986 ± 0.049 


Civ 


7a 


3.6+ 1.1 


171 ±9 


13 ±2 


0.892 ±0.126 


Civ 


7 


6.4 ±0.8 


232 ±6 


13 ±2 


0.892 ±0.126 



The scattered light component in this fit is 0.05 ± 0.01. 
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Table 10. Intrinsic Absorption Features in the STIS Spectrum of Mrk 509 



Ion 


Component 


N 


v sys 


b 


fc 




# 


(10 12 cm- 2 ) 


(km s _1 ) 


(kms" 1 ) 




Hi 


la 


11.6 ±4.2 


-428 ± 7 


17 ±2 


0.561 ± 0.052 


Hi 


1 


63.9 ± 14.5 


-380 ± 6 


32 ±5 


0.561 ±0.052 


Hi 


lb 


11.9 ±9.2 


-359 ± 6 


14 ±6 


0.905 ± 0.052 


Hi 


2 


459.2 ± 276.4 


-314 ±7 


24 ±6 


0.922 ± 0.019 


Hi 


2b 


86.8 ±73.1 


-283 ± 16 


27 ±12 


0.850 ±0.198 


Hi 


3 


87.5 ±41.8 


-237 ± 8 


34 ±4 


0.861 ±0.101 


Hi 


4a 


250.0 ± 98.3 


-66 ±6 


30 ±3 


0.960 ± 0.010 


Hi 


4 


180.3 ± 38.0 


-16 + 9 


18 ±6 


0.979 ± 0.026 


Hi 


5 


351.0 ±75.2 


37 ±6 


30 ±7 


0.958 ± 0.049 


Hi 


6a 


442.3 ± 188.9 


137 ±9 


43 ±2 


0.998 ±0.051 


Hi 


6 


401.4 ±709.1 


137 ± 26 


43 ±2 


0.998 ±0.051 


Hi 


6b 


6.8 ± 34.6 


168 ± 152 


43 ±2 


0.998 ±0.051 


Hi 


7a 


75.5 ± 16.2 


187 ± 134 


19 ±94 


0.632 ± 0.277 


Hi 


7 


206.7 ± 28.6 


220 ±7 


34 ±4 


0.632 ± 0.277 


Nv 


1 


95.5 ±34.1 


-408 ± 6 


10 ±1 


0.562 ± 0.094 


Nv 


lb 


16.6 ± 29.8 


-375 ± 6 


21 ±12 


0.810 ±0.104 


Nv 


2 


238.9 ± 56.2 


-319±7 


12 ±2 


0.860 ± 0.072 


Nv 


2b 


91.7 ±79.6 


-291 ± 6 


16 ±5 


0.809 ± 0.263 


Nv 


3 


94.0 ± 34.3 


-246 ± 6 


18 ±3 


0.759 ± 0.205 


Nv 


4a 


89.8 ±21.6 


-57 ±6 


14 ±1 


0.919 ±0.133 


Nv 


4 


263.1 ±23.3 


-17±6 


17+ 1 


0.997 ± 0.074 


Nv 


5 


44.0 ±25.9 


34 ±6 


18 ±2 


0.999 ± 0.057 


Nv 


6a 


16.7 ±7.1 


79 ±7 


41 ±6 


0.997 ±0.117 


Nv 


6 


46.8 ± 6.8 


121 ±6 


12 ±1 


0.997 ±0.1 17 


Nv 


6b 


5.2 ±3.0 


147 ±9 


11 ± 3 


0.997 ±0.117 


Nv 


7a 


25.1 ± 16.9 


184 ±8 


16 ±3 


0.704 ±0.231 


Nv 


7 


22.0 ± 14.8 


222 ±8 


16 ±3 


0.704 ±0.231 


Si iv 


4 


3.5 ±1.1 


-18 ±5 


6± 1 


0.946 ±0.197 


Civ 


1 


71.6 ± 19.6 


-406 ± 6 


12 ± 1 


0.560 ± 0.276 


Civ 


lb 


21.1 ±5.7 


-365 ± 5 


27 ±8 


1.000 ±0.110 


Civ 


2 


144.2 ± 24.2 


-320 ± 6 


14 ±1 


0.983 ± 0.078 


Civ 


2b 


137.9 ±23.3 


-292 ± 6 


30 ±5 


0.893 ± 0.090 


Civ 


3 


50.5 ± 27.8 


-243 ± 7 


25 ±2 


0.963 ±0.031 


Civ 


4a 


67.9 ± 8.5 


-60 ±6 


13 ±1 


1.000 ±0.069 


Civ 


4 


233.7 ± 17.9 


-18 ±5 


17 ±1 


0.960 ± 0.037 


Civ 


5 


35.8 ±45.6 


36 ±6 


20 ±1 


0.916 ±0.098 


Civ 


6a 


0.1 ±0.7 


88 ±68 


12 ±56 


1.000 ±0.045 


Civ 


6 


16.1 ± 18.1 


124 ±1 


13 ±1 


1.000 ±0.045 


Civ 


6b 


12.2 ±7.1 


120 ±1 


24 ±1 


1.000 ±0.045 


Civ 


7a 


1.1 ±1.2 


171 ±6 


6±3 


0.560 ± 0.276 


Civ 


7 


6.1 ±3.0 


214 ±2 


6±3 


0.560 ± 0.276 



The scattered light component in this fit is 0.01 ± 0.03. 
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Table 11. Intrinsic Absorption Features in the 1999 FUSE Spectrum of Mrk 509 



Ion 


Component 


N 


v sys 


b 


fc 




# 


(10 12 cm" 2 ) 


(kms- 1 ) 


(km s" 1 ) 




Hi 


la 


24.4 ± 9.9 


-433 ± 33 


32 ± 14 


0.991 ± 0.004 


Hi 


1 


363.3 ±9.1 


-411 ±7 


31 ±3 


0.991 ± 0.004 


Hi 


lb 


0.3 + 2.8 


-383 ±134 


15 ±99 


0.991 ± 0.004 


Hi 


2 


549.4 ±61.9 


-323 ± 6 


19 ±3 


0.978 ± 0.036 


Hi 


2b 


58.1 ±55.3 


-292 ± 5 


12 ±6 


0.826 ± 0.246 


Hi 


3 


681.5 ±104.9 


-258 ± 8 


31 ±3 


0.869 ± 0.037 


Hi 


4a 


1924.1 ±154.4 


-35 ±8 


43 ±2 


0.939 ± 0.029 


Hi 


4 


2594.9 ± 1146.8 


-29 ±7 


23 ±3 


0.902 ± 0.088 


Hi 


5 


359.5 ±45.2 


27 ±9 


32 ±9 


0.999 ±0.116 


Hi 


6a 


114.3 ±92.8 


75 ± 14 


32 ±1 


0.902 ± 0.020 


Hi 


6 


520.3 ± 93.7 


115 ± 6 


32 ± 1 


0.902 ± 0.020 


Hi 


6b 


168.4 ± 1.7 


129 ± 13 


32 ± 1 


0.902 ± 0.020 


Hi 


7a 


154.4 ±73.9 


174 ±8 


19 ±4 


0.905 ±0.154 


Hi 


7 


227.8 ± 118.6 


211 ± 8 


24 ±7 


0.969 ± 0.098 


Ovi 


la 


0.7 ± 5.5 


-427 ± 33 


9± 11 


0.697 ± 0.072 


Ovi 


1 


215.0 ±47.2 


-388 ± 8 


22 ±2 


0.697 ± 0.072 


Ovi 


lb 


154.9 ± 39.0 


-360 ±5 


22 ±2 


0.697 ± 0.072 


Ovi 


2 


566.2 ± 118.0 


-318 + 5 


19 ±3 


0.926 ± 0.076 


Ovi 


2b 


1248.1 ±395.5 


-290 ±5 


17 ±2 


0.919 ±0.037 


Ovi 


3 


675.2 ± 157.1 


-247 ± 5 


28 ±3 


0.814 ±0.048 


Ovi 


4a 


804.5 ± 387.2 


-65 ±5 


15 ±2 


0.714 ±0.049 


Ovi 


4 


8797.0 ±4120.3 


-9 ±5 


28 ±3 


1.000 ±0.008 


Ovi 


5 


683.5 ± 429.3 


47 ±5 


24 ±22 


1.000 ±0.032 


Ovi 


6a 


518.8 ±398.5 


87 ±5 


15 ±21 


0.949 ±0.160 


Ovi 


6 


3436.1 ±6985.0 


120 ±6 


12 ±7 


1.000 ±0.014 


Ovi 


6b 


104.5 ± 234.6 


146 ±5 


12 ±7 


1.000 ±0.014 


Ovi 


7a 


660.2 ± 176.7 


179 ±41 


22 ±78 


0.834 ±0.191 


Ovi 


7 


667.4 ± 1105.3 


213 ± 17 


26 ±8 


1.000 ±0.023 



The scattered light component in this fit is 0.04 ± 0.02. 

Table 12. Intrinsic Absorption Features in the 2000 FUSE Spectrum of Mrk 509 



Ion 


Component 


N 


v sys 


b 


fc 




# 


(10 12 cm" 2 ) 


(km s _1 ) 


(km s" 1 ) 




Hi 


la 


31.1 ±41.1 


-425 ± 8 


29 ±25 


1.000 ±0.071 


Hi 


1 


302.5 ±41.9 


-402 ± 5 


28 ±2 


1.000 ±0.071 


Hi 


lb 


0.3 ± 2.8 


-371 ±118 


15 ±19 


1.000 ±0.071 


Hi 


2 


563.0 ±167.1 


-318 + 5 


24 ±5 


0.970 + 0.078 


Hi 


2b 


75.9 ±69.2 


-281 ±8 


16 ±7 


0.890 + 0.072 


Hi 


3 


824.6 ±5 1.5 


-250 + 5 


28 ±1 


0.793 ±0.061 


Hi 


4a 


3493.7 ± 12430.4 


-35 + 54 


34 ±21 


0.904 + 0.093 


Hi 


4 


2506.3 ± 9025.3 


-4 ±48 


27 ±66 


0.793 ±0.159 


Hi 


5 


327.8 ±37.1 


47 ±53 


25 ±20 


1.000 ±0.027 


Hi 


6a 


175.9 ± 275.9 


89 + 107 


31 ±30 


1.000 ±0.027 


Hi 


6 


655.7 ±1019.0 


123 ±12 


31 ±30 


1.000 ±0.027 


Hi 


6b 


138.0 ±358.2 


166 + 138 


31 ±30 


1.000 ±0.027 


Hi 


7a 


365.8 ± 282.3 


185 ±6 


17 ±44 


0.565 + 0.394 


Hi 


7 


243.4 ±81.4 


217 + 26 


25 + 162 


1.000 ±0.027 


Ovi 


la 


16.5 ± 17.6 


-427 + 5 


9±5 


0.697 ±0.271 


Ovi 


1 


154.9 ±113.5 


-396 + 6 


17 ±7 


0.697 ±0.271 


Ovi 


lb 


210.5 ± 123.3 


-363 + 15 


17 ±7 


0.697 ±0.271 


Ovi 


2 


968.3 ±187.2 


-312 + 5 


18 + 4 


0.926 + 0.025 


Ovi 


2b 


125.6 ±68.6 


-281 ±7 


15 ±6 


0.919 ±0.103 


Ovi 


3 


917.3 ±181.2 


-256 + 5 


39 ±2 


0.814 ±0.065 


Ovi 


4a 


640.6 ± 68.3 


-56 ±5 


23 ±2 


0.924 + 0.008 


Ovi 


4 


1293.2 ±735.3 


-9 ±6 


20 ±1 


1.000 ±0.015 


Ovi 


5 


569.2 ± 275.9 


36 ± 10 


20 ±1 


1.000 ±0.052 


Ovi 


6a 


324.8 ± 179.7 


64 ±13 


20 ± 1 


0.949 + 0.082 


Ovi 


6 


565.4 ± 46.2 


104 ±5 


20 ±1 


1.000 ±0.027 


Ovi 


6b 


362.4 ± 59.3 


137 ±6 


20 ±1 


1.000 ±0.027 


Ovi 


7a 


403.8 ± 187.2 


185 ±6 


20 ±1 


0.834 ±0.160 


Ovi 


7 


640.3 ± 217.3 


213 ±5 


25 ±1 


1.000 ±0.019 



The scattered light component in this fit is 0.05 ± 0.01. 



